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Elements of the analytic structure of anomalous scaling and intermittency in fully developed hydro- 
dynamic turbulence are described. We focus here on the structure functions of velocity differences 
that satisfy inertial range scaling laws S n (R) ~ R l ’ n , and the correlation of energy dissipation 
K ee (R) ~ R-^. The goal is to understand the exponents and fi from first principles. In paper 
II of this series it was shown that the existence of an ultraviolet scale (the dissipation scale rj) is 
associated with a spectrum of anomalous exponents that characterize the ultraviolet divergences of 
correlations of gradient fields. The leading scaling exponent in this family was denoted A. The 
exact resummation of ladder diagrams resulted in the calculation of A which satisfies the scaling 
relation A = 2 — ( 2 - In this paper we continue our analysis and show that nonperturbative effects 
may introduce multiscaling (i.e. £„ not linear in n) with the renormalization scale being the infrared 
outer scale of turbulence L. It is shown that deviations from K41 scaling of S n (R) (£„ ^ n/ 3) must 
appear if the correlation of dissipation is mixing (i.e. p > 0). We derive an exact scaling relation 
fi = 2 £ 2 — £ 4 . We present analytic expressions for for all n and discuss their relation to experi¬ 
mental data. One surprising prediction is that the time decay constant T n {R) oc R Zn of S n (R) scales 
independently of n: the dynamic scaling exponent z„ is the same for all n-order quantities, z„ = 


1. INTRODUCTION 

This paper is the last in a series of four papers |]],|[[|. 
The aim of these four papers is to lay out the analytic ba¬ 
sis for the description of scaling properties in fully devel¬ 
oped hydrodynamic turbulence. We are concerned here 
with the statistical properties of turbulence in terms of 
averages over the fields of the fluid. The fundamental 
field in hydrodynamics is the fluid’s Eulerian velocity, de¬ 
noted as u(r, t) where r is a point in d-dimensional space 
(usually d = 2 or 3) and t is the time. The statistical 
quantities that have attracted decades of experimental 
and theoretical attention (see for example ||,||,|j,[7]]|] ) are 
the structure functions of velocity differences, denoted as 
S n (R) 

S n {R) = (|u(r + R,f) — u(r,f)| n ) (1.1) 

where (...) stands for a suitably defined ensemble aver¬ 
age. It has been asserted for a long time that the struc¬ 
ture functions depend on R as power laws: 

S n (R) oc R Cn , (1.2) 

when R is inside the so called “inertial range”, i.e. 
77 <C R <SC L. Here 77 and L are respectively the inner 
(viscous) and outer (energy containing) scale of turbu¬ 
lence. One of the major questions in fundamental tur¬ 
bulence research is whether the scaling exponents C, n are 


correctly predicted by the classical Kolmogorov 41 the¬ 
ory Q (known as K41) in which £„ = n/ 3, or whether 
these exponents deviate from n /3 as has been indicated 
by experiments. In particular we want to know whether 
the exponents may manifest the phenomenon of “mul- 
tiscaling” with £ n being a nonlinear function of n. 

Experimental research has not confined itself to the 
measurement of the structure functions of velocity differ¬ 
ences. Gradient fields have featured as well. For example 
the correlation function of the energy dissipation held has 
been studied extensively. The dissipation held e(r ,t) is 
dehned as 

e(M) = ^[V a ^(r,f) + Vf 3 U a (r,t)] 2 (1.3) 

with v being the kinematic viscosity. The correlation 
function of the dissipation held K ee (R) is 

Kee(R) = (e(r + R, t)e( r,t)) , (1.4) 

where e(r ,t) = e(r,f) — e. Here and below e = (e). Ex¬ 
periments appeared to show that K ee (R) decays 

according to a power law, 

K ee (R) ~ R~^ , 77<R<L (1.5) 
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with /r having a numerical value of 0.2 — 0.3. The an¬ 
alytic derivation of this law from the equations of huid 
mechanics and the calculation of the numerical value of 


the scaling exponent /.i have been among the elusive goals 
of theoretical research. 

In paper 0 of this series (subtitled “Line-Resummed 
Diagrammatic Perturbation Theory”) we reviewed the 
literature ( |||n],[yj with the aim of introducing the stu¬ 
dent to the available techniques. That paper did not 
present any new results. In paper I of this series || 
(subtitled “The Ball of Locality and Normal Scaling”) 
we dealt with the perturbative theory of the correlation, 
response and structure functions of hydrodynamic turbu¬ 
lence. The main result of that paper (and see also 1121) 
was that after appropriate resummations and renormal¬ 
izations the perturbative theory for these quantities is 
finite order by order. All the integrals appearing in the 
theory are convergent both in the ultraviolet and in the 
infrared limits. This means that there is no perturbative 
mechanism to introduce a length scale into the theory of 
the structure functions. In turn this result indicated that 
as far as the perturbative theory is concerned there is no 
mechanism to shift the exponents away from their K41 
values. Of course, nonperturbative effects may furnish 
such a mechanism and therefore in paper II || (subti¬ 
tled “A Ladder to Anomalous Scaling”) we turned to the 
analysis of nonperturbative effects. It was shown there 
that the renormalized perturbation theory for correlation 
functions that include velocity derivatives (to second or 
higher power) exhibit in their perturbation expansion a 
logarithmic dependence on the viscous scale 77 ||l3| , |l4| . 
In this way the inner scale of turbulence appears explic¬ 
itly in the analytic theory. The perturbative series could 
be resummed to obtain integrodifferential equations for 
some many-point objects of the theory. These equations 
have also non-perturbative scale-invariant solutions that 
may be represented as power laws of 77 to some exponents 
A. We argued p5t] that if A < A c where 


A c = 2 - C 2 


( 1 . 6 ) 


(a situation referred to as the “subcritical scenario”), 
then K41 scaling is asymptotically exact for infinite Re. 
In this case K ee (R) ~ e 2 (? 7 / 7 ?) 2 ( 4 / 3_A ), and the expo¬ 
nent /r is identified with 2(4/3 —A). The renormalization 
length is then the inner length 77 . 

In paper II it was shown that the exponen t A takes 
on exactly the critical value A = A c given by (|l. 6 | ), and 
the subcritical scenario is not realized. In the present 
paper it will be shown that the correlation function of 
the energy dissipation field has the dependence: 

K a (R) ~ e 2 (L/R) 2 C 2 " C 4 (I?/77) 2 ( A - A ^. (1.7) 


As a result of the criticality of A the inner scale 77 disap¬ 
pears from the correlation K ee (R) which finally takes on 
the form 


K ee (R) ~ e 2 (L/R) 2C2 " C4 . (1.8) 

Note that this implies that for K41 scaling (in which 
2 C 2 = C 4 )) Kee(R) does not decay as a function of R (i.e. 


the correlation is not mixing). This is hardly physical 
for a random field. On the other hand if there exists 
anomalous scaling in the sense that £4 < 2 ( 2 , then mixing 
is regained. Thus we will argue that the critical scenario 
A = A c goes hand in hand with anomalous scaling if 
K ee (R) is mixing, and then /i is identified with ((4 — 2 ( 2 , 

M = C4 - 2C 2 ■ (1.9) 


This is one of the exact scaling relations that are derived 
below. 

The main purpose of the present paper is to eluci¬ 
date the possible mechanisms for anomalous (non-K41) 
scaling of the structure functions in the theory of turbu¬ 
lence in fluids whose dynamics is described by the Navier- 
Stokes equations. We show that the analysis of the statis¬ 
tical theory of Navier Stokes dynamics leads to the possi¬ 
bility of multiscaling with nonlinear dependence of on 
n. For such a scenario we offer an evaluation of the scal¬ 
ing exponents C, n from first principles. We will show that 
the mechanisms for anomalous scaling are wholly non¬ 
perturbative. Order by order considerations always lead 
to K41 scaling. We show that a useful analysis of anoma¬ 
lous scaling can be developed on the basis of the so-called 
“balance equations” which are non-pertur bative. In sec¬ 
tion 2 we derive the balance equations ( 2.24 , 2.251 ) by 
writing the equations of motion for the structure func¬ 
tion S n (R ) and of related quantities. These equations, 
in the stationary state, exhibit a balance between a con¬ 
vective (interaction) term which is denoted as D n (R) and 
a dissipative term, denoted by J n (R), which is the cross¬ 
correlation between the energy dissipation e and (n — 2) 
velocity differences across a scale R. Both terms are ex¬ 
pressed as many-point correlation functions that depend 
on many coordinates, but some of these coordinates are 
the same. This last fact leads to the main strategy of 
this paper, which is to understand the fusion rules which 
describe the scaling of many-point correlation functions 
when some of their coordinates fuse together. When this 
happens the distance between the coordinates crosses the 
dissipative scale, and at that very moment the ultraviolet 
divergences are picked up. Since we learned in paper II 
how to compute these exactly, we have an important part 
of the fusion rules at hand. Following this strategy we 
succeed to compute exactly the dissipative term J n {R ) 
and the correlation function K tt {R). 

The fusion rules are developed in sections 3,4 and 5. 
In section 3 we show that the exponent A of paper II 
is indeed the relevant exponent for describing pair coa¬ 
lescence. In section 4 we discuss four-point correlations 
and their 2-point fusion rules. This leads to the calcula¬ 
tion of the correlation K ee (R) and the dissipative term 
J 4 . In section 5 we discuss the fusion rules of 7i-point 
correlation functions, leading to the exact evaluation of 
J n (R). A consequence of this calculation is the deriva¬ 
tion of the dynamical scaling exponents z n describing the 
characteristic decay time r n (R) oc R Zn of 77-point, 77 -time 
correlation functions of BL-velocity differences. We ar- 
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gue at the end of section 5 that all the exponents z n are 
independent of n and that they equal £ 2 - 

In section 6 we turn to the evaluation of the interaction 
term D n {R). Unfortunately we have not yet succeeded to 
evaluate D n exactly. Accordingly we describe in sections 
6 and 7 various possible evaluations of D n and the result¬ 
ing implications on £ n . It is shown in section 6 that the 
naive evaluation of D n leads to the scaling exponents of 
the /3-model. We give arguments however leading to the 
conjecture that the naive evaluation of D n is an over¬ 
estimate. We propose that there may exist a delicate 
cancellation in this evaluation, and the next order eval¬ 
uation leads unequivocally to multiscaling. In order to 
compute the scaling exponents £ n we need to know D n 
very precisely, coefficients and all. 

Not having an exact theory for D n we resort in sec¬ 
tion 7 to modeling. Guided by some diagrammatics and 
some common sense we conjecture an analytic form for 
D n which conforms with our expectations about eddy- 
viscosity. This part of the paper is not rigorous, and 
the results of section 7 should be considered therefore 
as tentative. The advantage of the approximation based 
on eddy viscosity is that we can offer analytic estimates 
of all the exponents ( n . We discuss these estimates and 
compare them with experimental findings in Table 1 of 
section 7B2. Section 8 is a summary and some discussion 
of the road ahead. 

Throughout the text we may refer to equations appear¬ 
ing in papers 0, I or II. When we do so we denote them 
as Eqs. (0-n), (I-n), (II-p) etc. 

2. DERIVATION OF THE BALANCE 
EQUATIONS 

Some of the most important nonperturbative con¬ 
straints on the statistical theory of turbulence are the 
balance equations (also known as the moment equations) 
for the structure functions. These equations are derived 
in this section and used in later sections to deduce the 
scenarios for multiscaling in the structure functions. Be¬ 
fore we derive the equations we introduce the statistical 
objects that appear naturally in the discussion. 

A. Structure functions and related quantities 

In this subsection we define quantities that are related 
to the structure functions S n (R). It is common in exper¬ 
iments to measure only the longitudinal structure func¬ 
tions [||. For theoretical treatment these quantities are 
not necessarily the most convenient since they are not 
invariant to rotations. It is useful therefore to intro¬ 
duce some related objects that have simple transforma¬ 
tion properties under rotations and inversions. The first 
one is the scalar quantity which is appropriate for even 


orders of S n . To keep in mind its scalar nature we will 

O 

denote it as S 2 m (R) and define it in Eulerian terms as 

S 2 m {R) = (|5u(r 0 |R,f)| 2m ) , R = r-r 0) ( 2 . 1 ) 

where 

<5u(r 0 |R, t) = u(r 0 + R, t) — u(r 0 , t) ( 2 . 2 ) 

is a simultaneous Eulerian velocity difference. The quan- 

O 

tity S 2 m (R) is analytic. For odd order structure func¬ 
tions we introduce a vector object S 2 m +i(R) according 
to 

^2771+1. C^-) — (5u a (r 0 |R,t)|du(ro|R, t)\ 2m ) . (2.3) 

Here and below we will use Greek indices to indicate vec¬ 
tor and tensor components, and Roman indices to indi¬ 
cate the order of the quantity. The placement of indices 
as subscripts or superscripts has no meaning, and is cho¬ 
sen for convenience. 

For isotropic turbulence the vector S% m+1 (R-) can only 
be oriented along R. This allows us the introduction of 
a scalar quantity S 2m +i(R) which depends on the mag¬ 
nitude of R: 

S2m+ i(R) = ^S 2m+ i(R) ■ (2.4) 

Lastly we will need also the tensor objects S^^R): 

^ 2 m+ 2 ( R ) = (<5w Q (r 0 |R,t)5w^(ro|R,f)|<5u(r 0 jR,f)| 2m ) . 

(2.5) 

Note that the objects introduced in (2T|, |2.3| . fHi| ) involve 
an arbitrary number of velocity fields but only two spa¬ 
tial points separated by R. We refer to them loosely as 
” two-point” correlation functions of velocity differences. 
The theory below calls for the introduction of three-point 
functions as well. It is best to define them using the four- 
point functions 

T“ r f +2 (R 1 ,R 2 ,R) = (fe Q (ro|Ri,t)^(ro|R 2 ,t) (2.6) 
x|<5u(r 0 |R, t)\ 2m ^ , 

^ 2 m +3 (Ri, R 2 , R) = (fe Q (ro|Ri,t)<Mro|R 2 ,t) (2.7) 

x/>M 7 (r 0 |R,f)|(5u(r 0 |R,f)| 2m ^ , 

r 2 m+ 4 (R 1 >R 2; R) = (5u a (ro\Ri,t)6uf,(r 0 \B*,t) (2.8) 

X(3u 7 (r 0 |R, t) 5 it 5 (r 0 |R, f)|<ju(r 0 |R, t)| 2m ^ . 

We will see that the theory produces expressions invlov- 
ing these quantities with two of the arguments being 
identical (i.e. Ri = R 2 , etc.) This fact will lead to 
the study of fusion rules in later sections. 
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B. The Balance Equations 

In this subsection we derive the balance equation that 
relates structure funtions to correlation functions involv¬ 
ing the dissipation field. We start by deriving some equa¬ 
tions of motion for the quantities defined above. 


1. Equations of Motion 

Our starting point is the Navier-Stokes equations for 
an incompressible fluid: 


where (V r ) 7 = d/dr 1 . Introduce now the shorthand no¬ 
tation w 2 and w 2m in situations in which the arguments 
are (r 0 |r,f): 

|w(r 0 |r,f)| 2 = w 2 , |w(r 0 |r, t)\ 2m = w 2m . (2.16) 


When other arguments appear they will be displayed ex¬ 
plicitlyIn term of these quantities the structure func¬ 
tions (2.1)-(2.3) are written as 


S2m (R) = ( W 2m ) , (2.17) 

S 2m+1 (I?) = (wu> 2m ) . (2.18) 


^ + [u(r, t ) • V]u(r, t) - u\7 2 u(r, t) - V p( r, t) 

= f(r, t ), V • u = 0 . (2.9) 

For simplicity we will choose the forcing such that V • 
f = 0. This equation can be rewritten in terms of the 
Belinicher-L’vov velocities v(jo,r,t) defined by Eq. (I- 
2.2) as follows: 

d V ^ r ’^ + [ w ( r °l r ’^ ' V]w(r 0 |r,- iA7 2 w(r 0 |r, f) 

- Vp(r 0 |r,t) = f(r 0 |r,t) , (2.10) 


where again R = r — ro and w without arguments means 
w(r 0 |r,f). Next observe that 


dw 2 


dt 


= 2mw 2(m “ 1) A 


dw 

~dt 


(2.19) 


Evalu ate now the scalar product obtained from Eq. 
( |2.15 ) by multiplying it on the left by 2mui 2 ( m ~ 1 )w. Us¬ 
ing Eq. ([2.19) this is written exactly as 


d S2m ( R ) 
dt 


+ D2m{R) — J2m(R) + Q2m(R ) j (2.20) 


where we also have the incompressibility condition V • 
w = 0. Here w is the BL velocity difference 

w ( r o|r, t) = v(r 0 |r,t) - v(r 0 |r 0 ,t) , (2.11) 

and p(r 0 |r, t), f(r 0 |r, t) are BL-transformed pressure and 
forcing (for more detail, see Section 3A of paper I). Ap¬ 
plying the transverse projector P this equation takes on 
the form 

+ P[w(r 0 |r,f) • V]w(r 0 |r,f) (2.12) 
- ^V 2 w(r 0 |r,f) = f(r 0 |r,t) . 


where we denoted 


D 2m (R)=2mJ dr 1 P a p(r 1 )-J^—(w 2< ' m 1} 


w 0 


ui 7 (r 0 |r + ri, t)in i a(ro|r + n,t) 

-w 7 ( r o|r 0 + ri,t)w; / 3(ro|ro + ri,t) (2.21) 


J 2m (l?) = 2 mv(w 2< ' rn l \ 


V 2 + V 2 

y r ' y r o 


( 2 . 22 ) 


Q 2 m{R) = ^rn(w 2{m 1] w a f a (r 0 \r,t) - f a (r 0 \r 0 ,t) 


(2.23) 


The application of P to any given vector field a(r) is 
non-local, and has the form: 


[Pa(r)] tt = / d,r'P a p(r — r , )ag(r / ), 


where P a p{ r — r') is 

P a p{ r - r') = 6 a p6(r - r') 


(2.13) 


(2.14) 


1 

4-7T 


S al 3 3 (U* - r'Jjrp -r'p) 


r / |3 


r f 15 


In deriving the equation for D 2m we used the incompress¬ 
ibility condition (which is not available in the Burgers’ 
equation) and performed changes of variables in the in¬ 
tegrals according to r' — r = iq and r' — r 0 = iq. In the 

O 

stationary state d 5 2m ( R)/dt = 0, and we can write the 
scalar balance equation 

D 2 m(R) = J 2 m(R) + Q 2 m(R ) • (2.24) 

Next we need to derive vectorial balance equations for 
S 2m -|-i. Repeating similar steps we end up with 


Next we consider Eq. (2.12) at two spatial points r and 
ro. Subtracting the two equations we get 

dw a (ro|r ,t) + J dr ^p a0 ( r _ r /) _ P a p(r 0 - r')] 

xw 7 (r 0 |r',t)(V t ./) 7 M; / 3 (ro|r , ,t) (2.15) 

= i/[V 2 + V 2 Ju> a (ro|r,f) + / a (r 0 |r,t) - /a(r 0 |r 0 ,t) , 


D2m+l(R) = J2m+l(R) + Q2m+1 (R) , (2.25) 


where 


D* m+1 {K) = 2rn J driP a/3 (ri)^-^ 


W 2 ( m d Ws yj a 


w 7 (r 0 |r + r 1 ,t)w f 3 (r 0 \r + r 1 ,t) 


4 






























-w 7 (r 0 |r 0 + r 1 ,t)w / 3 (r 0 \r 0 + r 1 ,t) 


+Jdr 1 P s/ 3 (ri) (w 2m [w 7 (r 0 |r + ri, t)w f 3 (r 0 |r + ri, t) 
-w; 7 (ro|ro+ ri,f)w /3 (ro|r 0 +ri,t)]^ , (2.26) 

J 2m+i( R ) = 2rnv(w 2( ' m - 1) w 5 w a [Vl + V 2 o ]w a ^ 

+v(w 2m [V 2 r +V 2 ro ]w s ) , (2.27) 

Q 5 2 m +1 ( R ) = 2m(w 2( ' m ~ 1) w s w a [fa (r 0 |r, t) - f a (r 0 |r 0 ,t)]j 
+(w 2m [f8(r 0 \r,t) - fs{r Q \r 0 ,t)]j . (2.28) 


For isotropic turbulence all these vector quantities have 
ob viou s scalar counterparts, similar to the one introduced 
in (2.4). For example, 


^2 Q m+ l( R ) = -jjr D 2in+l(R) • , (2.29) 

^2m-\-\ ( R ) = ^J2m+l(R) • (2.30) 


With thes e definitions we can rewrite the balance equa¬ 
tion (2.25|) in the form 


D2m+l(R) = J2m+l(R) + Q2m+l(R) ■ 


(2.31) 


Together with ( 2.24 ) we can finally write for n odd or 
even 


other hand the Burgers equation describes a compressible 
flow, and this fact introduces the coefficient (n + 1 )/n 
in the RHS of this expression. For the Navier-Stokes 
equations we can only use the expression for a rough 
evaluation of D 2m . Notwithstanding, it is exact for the 
case m = 1, where it takes the form 

D 2 {R) = ^-S%{ R) . (2.35) 


The deep reason for this simplification in the case to = 
1 is the conservation law (in the absence of viscosity) 
of the quadratic invariant of the BL equation of motion 
( M0fc , i-e. / drw 2 [r 0 |r,f]). D 2 (R) is simply the rate 

of change of w 2 [r 0 |R, t] and therefore must be a pure 
divergence. None of the higher powers of w[ro |R, t] forms 
an invariant, and as a consequence none of the higher 
orders D n (R) can be written as a pure divergence. On 
the other hand we see from Eq. (2.21) that with the full 
kernel we have two terms that can be expressed in terms 
of the T tensors: 


D 2rn (R) = 2m J dr 1 P a p(r 1 ) (2.36) 

X [T^ +1 (R + n, R + iq, R) — T^ +1 (n, r,, R)] . 

This is the final form of D 2m (R). 

Following the same steps of analysis we find the final 
expression for D 2 m +i(R). 


D n (R)=J n (R) + Q n {R). (2.32) 

We discuss each term in the balance equation in the fol¬ 
lowing subsections. 


C. The Terms of the Balance Equation 

1. The Interaction Term 


R^2m+l (R) = 2m J dr 1 P a p(r 1 ) (2.37) 

X ^{«( R + rr,R + n,R) - T^n,n, R)} 

+ J driP 5 / 3 (ri)^-|r 2 / 3 J +2 (R + ri,R-(-ri,R) 

—T 2 ^ +2 (ri, ri ,R)} . 


To understand how to evaluate D 2m {R) we begin the 
discussion by neglecting the effect of the pressure term 
in the Navier-Stokes equation. This is equivalent to re¬ 
placing the kernel of the projection operator P a ff{r ) by a 
the unit operator <5 Q , ) g<5(r). The effect of this replacement 
of kernels is to change D 2m (R ) to 

D2m(R) = 2m ^w 2 ^ m ~ 1 ' ) w a w ry ^^-'^ . (2.33) 

This may be rewritten as 

D2m(R) = {w 1 w 2m ) = ^-^ m+1 (R) (2.34) 

(neglecting pressure) . 


The naive evaluation of every term in Eq. (2.36) is given 
by (2.34). However this evaluation is only acceptable 
under two assumptions: (i) that the integral converges, 
and (ii) that there are no cancellations between the terms 
with opposite signs in Eq. ( 2.3C| ) . The analysis of conver¬ 
gence requires understanding the asymptotic properties 
of the T correlators. Similar properties will determine the 
evaluation of J 2m as will be seen next. For that reason 
we will devote the next Sections to questions of asymp¬ 
totics and fusion rules for n-point correalation function 
of ve locity differences. We will learn that the integral in 
Eq. ( 2.36| ) indeed converges, but on the other hand can¬ 
cellations are not excluded and are the possible source of 
multiscaling in turbulence. In fact, such a cancellation is 
the only mechanism of multiscaling that we can identify. 


This expression is relevant for the 1-dimensional Burgers 
equation, since there is no pressure in that case. On the 
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2. The Dissipation Term 


Starting with (2.22) we first use the fact that for simul¬ 
taneous correlations we can use the Eulerian differences 
(5u instead of the BL-velocity differences. Secondly we 
use the symmetry with respect to exchange of r and ro 
to write the expression with twice the Laplacian with re¬ 
spect to ro only, instead of the sum of Laplacians. Lastly 
we use translational invariance to move one of the gradi¬ 
ents around. We find 


J 2 m{R) = -4mz/< ( |(Su(r 0 |R)| 2(m 11 


d 


drop 


Ua(r o) 


+2 (m - 1)( |5u(ro|R)| 2(m “ 2) fe ct (ro|R)<5u 7 (ro|R) 


drop 


r( r o) 


dr 0 p 


t( r o) 


(2.38) 


To relate this quantity to the tensors T defined in Eqs. 
(2.6-2.8) we write the derivative as the limit of velocity 
differences: 


■7T— ««(r 0 ) = lim 

dr 0 p dp^o 


u a { r 0 + dpep) - u a (ro) 
dp 


(2.39) 


where ep is the unit vector in the f3 direction. With this 
in mind the expression for J 2 m is 

1 


J 2 m{R) = —4 mu lim l0 

d,3^0 1 


T 2 m(dpep,dpep,R) 

+2(m-l)T^{dpep,dpep,R)] . (2.40) 


One should understand that a summation over a, /? and 
7 is implied as usual. Note that for dp <C rj the velocity 
field is expected to be smooth, and both T functions 
become proportional to t#g, cancelling this factor in the 
denominator. We will not have a dependence of the limit 
on the direction of the vector d = dpep. 

It can be seen that the evaluation of J 2 m{R) again 
requires elucidation of the asymptotic properties of the 
correlations functions including the rules of coalescence 
of groups of points. In paper I we proved the property 
of locality that means that the coalescence itself is reg¬ 
ular. However, we have here differential operators on 
coalescing points and we need to consider not only lead¬ 
ing behaviours but also next order properties. This is 
done in the section 5. 

The quantities T|' rn+1 (R) can be expressed in terms of 
the third and fifth rank T tensors, but we avoid display¬ 
ing the result since we do not use it explicitly. 

Lastly we need to discuss the forcing terms Q 2 m{R) 
and Q£ m+1 (R). This is the topic of the next subsection. 


3. The Forcing Term 


finds this statement believable can skip this subsection. 
To discuss the forcing terms it is useful to introduce a 
higher order Green’s function according to 


Q 


a(3 

2m—1 


,1 


(yo\x,x') 


I dw 1 d n ^Wg \ 
\ 6fp(r 0 \x') / 


(2.41) 


We remind the reader that x = (r, t) and that according 
to our convention we can omit the argument (ro[ir) in the 
function w. 

When the statistics of the r ando m forcing is Gaussian 
we can write Q 2 m(R) of Eq. ( 2.22 ) as 


Q2m(R)= I dr'dt'g^_ 1A (r Q \x,x') D a p(r 0 \r', r, t' - t) 


-Dgp(ro\r',r 0 ,t' - 1) 


(2.42) 


Here, as usual, R = r r 0 and D a p( r 0 |x, a;') is the cor¬ 
relator of f a{x) and fp(x). 

We are going to estimate Q 2m (R) by taking first a 
delta-correlated forcing, 

D a p{vo\x, x') = Dgp (r - r')8{t - t') . (2.43) 


The consideration of the effec t of f inite c orrel ation time 
is deferred to the end. Using ( 2.43| ) Eq. ( 2.42p simplifies 
to 


Q 2 m{R)= I dr'g^_ 1A (r o \r,r\0) D a0 { r'- r) 


- D a p(r' - r 0 ) 


(2.44) 


Because of the delta functions we have lost the time in¬ 
tegration and we have the zero-time Green’s function, 
which is not dressed by the interaction (and see paper I 
for an explicit de mons tration). Using the chain rule of 
differentiation in (|2.4l|) we get 




U = / w 2 ^ 


-i). 


Sw n 


+2(m - 1) ( w 2{m - 2) wgw. 


8fp(r 0 \x') 
6 w-/ 

8fp( r o\x') 


(2.45) 


Next we recall that at time t = 0 the unaveraged response 
Sw/5f is uncorrelated with the velocity field since any in¬ 
teraction involves a vertex with time integration between 
0 and t. Accordingly we can decouple the response in 
( p.45 ) and write 

= 0°/?( r o|r,r') S 2 (m- 1 ) (R) (2.46) 

+2(m-l)^ /3 (ro|r,r , )S , 2 ( L-i)( i? ) . 


In this subsection we show that the forcing term can be 
neglected in the inertial range of scales. The reader who 
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where g° is the bare Green’s function which was defined 
and computed in paper I sections 2,3. We proceed to 
evaluate the order of magnitude of Q 2 m by substituting 

































the last equation in ( [2.44 ). Suppressing vector indices we 
evaluate 


3. THE ANOMALOUS EXPONENT 
ASSOCIATED WITH PAIR COALESCENCE 


Q 2m (R) ~ m(2m - l)S 2(m _ 1) {R)Q 2 {R) , (2.47) 

where 


Q 2 (R) = / dr'(7 o (ro|r,r',0) D(r' - r) - D(r' - r 0 ) 


(2.48) 


We can estimate Q 2 by recalling that D{R) has a char¬ 
acteristic scale L (which is the outer scale of turbulence). 
For R L we expand D(R) in R: 

D(R) = D( 0) + D"R 2 /2 , D" ~ D(Q)/L 2 . (2.49) 


In its turn D( 0) is about the rate of energy input, and 
for stationary turbulence D( 0) ~ e, the rate of energy 
dissipation. Finally we estimate 

D{R) - D{ 0) ~ e(R/L) 2 . (2.50) 

Together with the evaluation of G° that was presented 
in paper I Section 3 B we evaluate Q 2 as 


Q 2 {R) ~ 


L 2 


dr' 


|R-r'| 2 

|R-r'| 3 



(2.51) 


The integral can be evaluated by introducing the ultra¬ 
violet cutoff at the viscous scale 77 : 


Q 2 


In 


(2.52) 


Substituting in (2.47) and omitting numerical factors we 
end up with 


Q 2 m(R) 


e 



2 

R‘2(m — 1) ( R ) 



(2.53) 


We are now in a position to comp are Q 2 m with the in¬ 
teraction term D 2 m given in ( 2.33| ). For K41 scaling and 
delta-correlated forcing we have 


Qim (R\\ fR\ 

5^ ~ UJ ln (?) 


(2.54) 


It should be noted that the logarithm in this equation is 
a direct result of our use of delta correlated forcing. This 
choice brought in the zero-time Green’s function with 
its inverse cubic R dependence, which led to the loga¬ 
rithm. Any realistic forcing would be integrated against 
the finite time response which is less singular than the 
zero time function. We thus assert at this point that a 
more detailed analysis should reveal that Q 2 m is negligi¬ 
ble compared to D 2 m in the core of the inertial interval. 
The power law ( R/L) 2 for the ratio of these quantities 
is appropriate for K41 scaling, and the exponent may be 
slightly different for non-K41 exponents. At any rate we 
are satisfied that for ]J<L the forcing term is negligible 
compared to the transfer term. 


In this section we begin the exploration of the ultravi¬ 
olet properties of n-point correlation functions with the 
aim of computing J n and D n of the previous section. 
To achieve this we examine first a second-order Green’s 
function of the type introduced in paper II and whose 
ultraviolet properties are characterized by an anomalous 
exponent A. We will prove that these ultraviolet prop¬ 
erties are shared by the second order structure function, 
and that A attains the critical value A = 2 — ( 2 - Al¬ 
though this statement was made already in paper II sec¬ 
tion 5.B, we discuss it here in full detail because of its 
crucial importance for the structure of the theory of tur¬ 
bulence, and in particular in determining also the ultravi¬ 
olet proper ties o f the quant ities appearing in the balance 
equations (2.24) and (2.25). We begin with the 4-point 
Green’s functions. The analysis of this section will al¬ 
low us already to eva luate the correlation function of the 
dissipation field (|l.4| ) and to derive the important scaling 
relation ( |l .9| ) . The strategy of this section is to first intro¬ 
duce the needed many-point Green’s functions, then to 
examine and classify their diagrammatic representation, 
and lastly to resum the diagrammatic series to obtain 
exact integral equations for these quantities. The exact 
equations will allow us to find non-perturbative solutions. 


A. Many-Point Green’s Functions 


1. Definitions 


Consider the nonlinear Green’s functions Gm,n which 
are the response of the product of m BL-velocity differ¬ 
ences to n perturbations. In particular 


£ 2 , i y (, r o\xi,x 2 ,x 3 ) 
£“,2 7 ( r o|a:i) * 2 , 2 : 3 ) 


S[w a (r 0 \xi)wf3(r 0 \x 2 )] \ 
6h 7 (x 3 ) / 

/ S 2 w a { r 0 |xi) 

\ Shp(ro\x2)Sh^ (x 3 ) 


(3.1) 

(3.2) 


(-> 06(378 

^2,2 


(r 0 |a;i, x 2 , £ 3 , * 4 ) 


/ 5 2 [wq(r 0 |a:i)w l 3(ro|^2)] 
\ 5h 1 {xz)5hs{x^) 


QzT / \r 0 \xi,X2,X3,X4,) 


_ / i5[w a (r 0 |.Ti)w /? (ro|a:2)w 7 (a:3)] \ 
5hs(xi) / 


(3.3) 

(3.4) 


Note that these are different from the nonlinear Greens’ 
function G 2 that was dealt with extensively in paper II, 


(Y 0 \xi,X2,Xz,Xi) 


i5ti; a (r 0 |xi) 5w 0 {tq\x 2 ) \ 
Sh 1 (xz) Shs{x 4 ) / 


(3.5) 
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One can see that the relation between these 4-point 
Green’s functions follows from the chain rule of differ¬ 
entiation and is 


, x 2,X 3 ,X4) = G 
+ G2 /3S ' 1 (r 0 \x 1 ,x 2 ,x 4: ,X3) + ^ 


%^\r 0 \xi,X2,X3,X4) 

u>a(xo\xi) 5 2 wp(r 0 \x 2 ) 

Sh- y (x3)Shs(x4 : ) 


+ 



/3(r 0 |a:2)^ 2 Wc t (ro|a;i) \ 

5h~ f (x 3 )6hs(x 4 ) / 


(3.6) 


In order to analyze the properties of the functions Gm,n 
we will explore their diagrammatic representation. 


2. Diagrammatic Representation 

In section 3 A of paper II it was explained in detail 
how to produce the diagrammatic representation of G' 2 . 
Very similar steps lead to the diagra ms f or C/ 2 , 2 , which are 
shown in Fig. [l. Of course, due to ([ht]) all the diagrams 
of G 2 appear also in the representation of 1 / 2 , 2 - These 
common diagrams all have two principal paths made of 
Green’s functions as explained in paper II. In G 2 2 ther e 
are contributions coming from the last two terms in ( |l. 6 | ) . 
To understand the structure of these new diagram note 
that every diagram representing 5 2 w\/ 5h 3 8h^ has a prin¬ 
cipal path of Green’s functions that starts with coordi¬ 
nate x\ or x 2 which splits at some point into two principal 
paths made of Green’s functions ending up with the co¬ 
ordinates £ 3 , X 4 . The field w 2 itself is another tree of the 
type appearing in Fig.3 of paper II. Upon multiplying 
these trees and averaging over the Gaussian ensemble of 
forces we get diagrams that are ready to be resummed. 
One set of diagrams that appear are those that have a 
bridge made of a single Green’s function connecting the 
“left” part of the diagram (coordinates ri, r 2 to the “right 
’’part of the diagram (coordinates r 3 ,r 4 ). By necessity 
this Green’s function belongs to the principal path of the 
diagram which is made of Green’s functions. We note 
that a bridge made of a single propagator may appear 
only once. A second singly propagator bridge makes the 
diagram “one-eddy reducible”, and such diagrams have 
already been line-resummed; see paper 0 for more de¬ 
tails. Every decoration to the left of the bridge may be 
resummed into a dressed vertex on the left, and every 
decoration on the right of the bridge may be resummed 
into a dressed vertex on the right. The fully resummed 
series of weakly linked (via one-propagator bridge) dia¬ 
grams which we denote as Q 22 contains exactly the two 
diagrams that are shown in Fig. [I], panels b-d. 

The reader’s attention should be drawn to the weakly 
linked diagrams which have appeared here for the first 
time. It will turn out that these diagrams play a very 
important role in the mechanism that we propose for 
anomalous scaling in turbulence. For that reason we 
pause for a moment to discuss the topological structure 
of these diagrams in more detail. By construction the 


bridge between points (1,2) and (3,4) in the diagrams for 
^ 2(2 x 2l x 3: £ 4 ) can consist of a Green’s function 
but not of a correlator. In panel c of Fig. [l] the coordi¬ 
nates of the Green’s function are denoted by x a ,Xb . On 
the two sides of the ( a — b) bridge there exist three-point 
objects. These objects have an exact resummed form in 
terms of the dressed vertices A and B which were briefly 
introduced in paper I. Note that we have a freedom in 
the diagrammatic representation as to whether to include 
the bridge itself together with the object on the right 
or on the left. In the former case the resulting object 
on t he right is Gi,2(xo\x a , X3, £ 4 ) which was defined in 
( |3.2D . In the latter case, which is the convention that we 
choose, see panel b, the resulting object on the left side 
is G- 2 ,i( r o\xi, x 2 , Xb) which is defined in ( |3.l| ). In order to 
display the bridge (a — b) explicitly we introduce in panel 
c a new object denoted as V(r 0 \xi, x 2 , x a ). This object 
has two ’’entries” 1 and 2 starting with propagators and 
one “exit” denoted x a which ends with a vertex, panel d. 
We will use an empty small circle to denote the position 
of a vertex. This is done to distinguish a vertex from a 
propagator leg. 

Next, as explained in paper II, we need to identify two- 
eddy reducible and two-eddy irreducible diagrams. The 
first type are diagrams that can be split into a left and a 
right part by cutting two propagators. For G 2 these prop¬ 
agators can only be Green’s functions. This is the cross 
section denoted by (a) in Fig. ||. In the present case we 
can also have cross sections of type (b) with one Green’s 
function and one propagator. There cannot be cross sec¬ 
tions with two correlators because we always have at least 
one principal path of Green’s functions that connects the 
left part of the diagrams to its right part. Finally we 
can have a cross section of type (c) in which two Green’s 
functions appear in opposite orientations. 

In summary, the structure of the diagrammatic series 
can be described as follows: we get ladder diagrams that 
have alternating Green’s functions and double-correlators 
in any possible order, as long as there is at least one 
Green’s function between two rungs which carries the 
principal path. The ladder diagrams can have alternat¬ 
ing rungs of three types, depending on the type of prop¬ 
agators preceding the rung. The three types of rungs are 
sums of two-eddy irreducible diagrams that are denoted 
S( 13 ), U( 2 , 2 ) and E( 3 j i) respectively. Graphically they 
are represented by a dashed bar, empty bar and doubly 
dashed bar respectively. The first one has three straight 
tails and one wavy tail, the second one has two straight 
and two wavy tails and the last has three wavy and one 
straight tail. The series for are shown in Fig. |] 

and Fig. [|. Fig.2 shows the diagrams for E( 2i2 ). This 
series is composed of all the diagrams in Fig.10a of paper 
II (and diagrams 1, 2 and 5 in Fig. |] are examples of 
those) in addition to new diagrams like diagrams 3 and 4 
in Fig. |^. All the old diagrams have a horizontal princi¬ 
pal cross section that cuts through correlators only. The 
new diagrams contain the split in the principal path, and 
the principal cross section through correlators turns up 





or down by 90°. The series for the two other two-eddy 
irreducible mass operator £(i 3 ) and £( 31 ) are shown in 
Fig. @ The diagrams in these series have no principal 
cross section that cuts through correlators only. 


B. Calculation of A 

In order to evaluate A we establish a fundamental iden¬ 
tity which is 


3. Resummation of the strongly linked contributions 


5F a p(r 0 \xi,x 2 ) 
5D lS (x 3 , x 4 ) 


s->a/3 , y8 

^ 2,2 


(r 0 |xi, £ 2 , X 3 , X 4 ) , 


(3.8) 


As in the case of the ladder diagrams of G 2 we can 
resum also in the case of G 22 all the ladder diagrams into 
integral equations. The resummed ladders are shown in 
Fig. | The diagrammatic notation of G 22 is an empty 
circle with two wavy and two straight lines. It has three 
different types of contributions. First come the reducible 
contributions that we denote by G 2 ° 2 and are represented 
by the unlinked diagrams (1) and (2) of Fig. |[ The first 
two ladders (3) and (4) in Fig. |]are identical to the RHS 
of the resummed equation for G 2 . The next term on the 
RHS in Fig. |] are new. The resummation of the ladders 
is shown in Figs. |5| On the RHS of the diagrammatic 
equation we find a Green’s function which is defined in 
(3.4) as G 3 ,i and whose graphic notation has a crossed 
circ e with three wavy and one straight tails. This object 
is again resummed in terms of itself and C ? 2j2 as shown in 
Fig.5b. 

In conclusion we learn that the 4-point Green’s func¬ 
tions Gm,n satisfy equations that contains linear opera¬ 
tors acting on Gm,n and inhomogeneous terms that are 
products of G’s and the weakly linked contributions. To 
be sure, the linear operators are themselves functionals of 
Gm,m so that our equations are in fact nonlinear. If we 
expand the solutions of these inhomogeneous nonlinear 
equations around the inhomogeneous terms we generate 
back the initial diagrammatic expansion. However, now 
we can also explore the “ homogeneous ” solutions of these 
equation that are obtained by discarding the inhomo¬ 
geneous terms. Such solutions, if they exist, are mani¬ 
festly nonperturbative effects. We will need to show a 
posteriori that these nonperturbative solutions are much 
larger than the solutions that can be found from per¬ 
turbative analysis. Indeed, in paper II it was shown in 
detail that G 2 has such a nonperturbative homogeneous 
solution with the following property: when the first two 
coordinates r\ and r 2 are of the same order (say r) and 
much smaller than the last two coordinates r 3 and r 4 
(which are of the order of R ) then 


where the covariance D is the correlation of the pertur¬ 
bations, 

D 1 s{x 3 ,xa) = (h 1 {x 3 )hs{x 4 )) . (3.9) 

The identity is proved most easily using the path integral 
formulation as reviewed in paper I. In terms of the func¬ 
tional Z( 1, m) of Eq. (1-3.12) the second order Green’s 
function is 

G 2 P 2 \r 3 \x i,x 2 ,x 3 ,x 4 ) 

= - (w a {ro\xi)wf}(Y Vl \x 2 )p 1 (x 3 )p 8 (x 4 )) . (3.10) 

On the other hand we see from Eq. (1-3.14) that the 
derivative with respect to D a 0 brings down ip a P0 '■ 

SD^hy) = ip -’ ( ' X}pM ■ (3 U) 

This means that the functional derivative of (wawp) with 
respect to D 7 s is precisely the RHS of ( [l.lCj ). This is the 
proof of the identity. 

Rewrite now the identity in the form 

5F af 3 (r 0 |xi,x 2 ) 

dx 3 dx 4 G 2 ^ r 6 (r 0 \x 1 ,x 2 ,x 3 ,X 4 ) 6 D j 5 (x 3 ,X 4 ) . (3.12) 

In this form this is a relation of the response SF in the 
velocity correlator F to a variation in the correlator of 
the random forcing SD. It is clear now that if the random 
forcing is limited to scales r 3 ,r 4 ri, r 2 , the existence of 
flux equilibrium with a scaling solution for F a p( ro|xi, x 2 ) 
means that the variation 5F a p(ro\x\,x 2 ) must be propor¬ 
tional to F a p( ro|xi,x 2 ) itself. We can have a change in 
the amplitude but not in the functional form: 

$F a p( r 0 |xi,x 2 ) oc F a / 3 (r 0 |xi,x 2 ) 

for n,r 2 -C r 3 ,r 4 . (3.13) 


V(iG 2 (r, r, R, R) ~ r~ A for R > r . (3.7) 

Since the same terms appear in the equation for C / 2j2 we 
know that such a divergence (with r —> 0 ) must appear 
also in all G n ,m- (This of course hinges on the assump¬ 
tion that the new terms in Fig. do not contribute an 
exact cancellation; in view of their different nature we 
judge this possibility unlikely). The new terms may have 
a stronger or a weaker divergence, and we turn now there¬ 
fore to an exact calculation of the value of A. 


To understand what are the requirements of the variation 
SPrys(x 3 , X4 ) tha t guarantee the validity of the universal 
behaviour ( 3.13 ) let us write again the Wyld equation 
(II—2.9). We will use economic notation, such that the 
vector index carries implicitly also the position coordi¬ 
nate. Repeated indices must be summed upon and the 
convention is that this sum also requires integration over 
the intermediate position coordinates. This convention 
allows us to write the Wyld equation as 


F a 0 — G(x'y[D'yS T ^ 7 s\Gf05 


(3.14) 
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The variation SD^s causes a variation in F. In our con- C. Three Point Objects and the Weakly Linked 

vention inhomogeneous Contributions 


F a p + SF a p — G ai [D 7 s + SDjs + . (3.15) 

Next note that $ 7 < 5 (ro|a;Q, Xb) can be exactly expressed 
as a second derivative of the 4-point correlation function 
F±(x a , x a , Xb, Xb), cf. Eq. (1-4.5). Consequently 

$ 7 5( r o \x a ,Xb) ~ r c a l ~ 2 . (3.16) 

We know that £4 is expected to be considerably smaller 
than 2. (The K41 estimate is £4 = 4/3, whereas ex¬ 
perimentally one finds (4 ~ 1.2). Thus d) 7 ^(ro|x a , Xb) is 
growing when the coordinates become smaller, whereas 
1 5Djs is restricted to the large scales and is decaying for 
smaller coordinates. We expect therefore that 6 F will 
be proportional to F and the constant of proportionality 
is determined by the boundary conditions at large scales 
where both D and SD are not negligible. 

For future reference we should note at this point that 
the proportionality of 6 F a P and F a P is not restricted 
only to their scaling exponents. In fact the two quantities 
have the same tensor structure. In other words they are 
the same function of ri — ro and r 2 — ro up to constants 
that are independent of the tensor indices. To complete 
the argument choose now t\ = f 2 - Next notice the fact 
that 


Vi-V 2 i^ 9 (ro|ri,r 2 )ocr& 2 . (3.17) 


Now restrict 5D 1 s{ £ 3 , 24 ) to , rq ~ R n,r 2 . Apply¬ 

ing the operator Vi • V 2 to ( |3. 12| ) we conclude that 

Vi • V 2 t7£ 2 7<5 (r 0 |ri, r 2 , X3, £4) oc r[l~ 2 . (3.18) 


Note that the RHS is a function of 7*12 only; the reason 
for this is that under the derivatives the dependence on 
ro disappears. Accordingly the restriction of validity of 
this result is not necessarily n, r 2 Cl? but just ri 2 -C R. 
Comparing with (3.7) we reach the central result of this 
section: 


A = 2 - C 2 • (3.19) 

Finally we can argue that the application of the oper¬ 
ator Vi ■ V 2 to Qs'i will give rise to the same exponent 
A as in ( |3.18| ). To see this apply the operator to the two 
equations in Fig.5. Suppose that the divergence associ¬ 
ated with C/ 3,1 is stronger than A. This will immediately 
force the divergence of the LHS of Fig.5a to be stronger 
than A, in contradiction with our exact result. In the 
same manner the divergence of Vi ■ V 2 C/ 3 ,i cannot be 
weaker than A due to the equation in Fig.5b. We thus 
conclude that also 

Vi • V 2 ^ 3)1 75 (r 0 1r 1 , r 2 , £ 3 , £ 4 ) a r^ 2 ~ 2 . (3.20) 


Equation ( 3.19| ) was derived by neglecting the inhomo¬ 
geneous part of the equation for Q. In order to be sure 
that the inhomogeneous part is not important we need to 
evaluate now the wea kly lin ked c ontributions and com¬ 
pare them with HI and ( |3.20| ). We are not going to 
perform the evaluation with the same care as we did in 
the computation of A. All that we need is to show that 
these terms are negligibly small. We will show this by 
assuming that the diagrams have the property of rigid¬ 
ity which was demonstrated in paper II order by order. 
This will be sufficient since we will be able to show that 
there exists a large gap in the value of the scaling expo¬ 
nents of the homogeneous and the inhomogeneous terms. 
Such a large gap is not expected to be swamped by non- 
perturbative effects. 

To prepare for this evaluation we need first the dia¬ 
grammatic representation of three-point objects. These 
are the 3-point correlator F 3 , and the Green’s functions 
C/ 2 p and C/i, 2 . All the three point objects can be repre¬ 
sented with the help of the three types of dressed vertices 
which were denoted in paper II as A, B and C respec¬ 
tively, see paper II Fig. 7. We remind the reader that 
vertex A is a junction of one straight and two wavy lines, 
vertex B is the junction of two straight and one wavy 
line, whereas vertex C is the junction of three straight 
lines. The diagrammatic representation of the triple cor¬ 
relation function Fj in terms of these vertices is shown 
in Fig. 7 of paper II, and is reproduced in more com¬ 
pact form in Fig. panels b,c. In the same manner the 
Green’s function Q 2 p is represented here in Fig. [j] panel 


c. 


Consider now the weakly linked diagrams for t/ 2j2 
which are resummed into the form shown in the diagram 
in Fig. |l] panel b. In the limit of r± ~ r 2 ~ r —> 0 and 
^3 ~ r 4 ~ R, rigidity means that the integral over rt 
contributes mostly in the region ~ R. We thus need 
to understand the r dependence of f? 2> i(ro|xi, x 2 , Xb) in 
the asymptotic situation ri,r 2 <C R. Looking at Fig. |l] 
panels c, d we see that due to rigidity the vertices A and 
B will contribute mostly in the region r a ~ r. In K41 
scaling the exponent of r is found by noticing that the 
correlator contributes r*' 2 , and the Green’s function to¬ 
gether with the integration over x a contribute r z . The 
Green’s function of the bridge a — b in Fig. [l] panel c 
contributes r/i? 4 , or r 1 . The vertex A itself in diagrams 
1 and 2 in Fig. [l] panel d gives r _1 , and in total we find 
r ( 2 +z f or these contributions. The same arguments give 
the same asymptotic behaviour of diagram 3 shown in 
Fig. 0 panel d. This behaviour should be compared with 
r ^ 2 for the homogeneous term, justifying the neglect of 
the inhomogeneous contribution. Note that small correc¬ 
tions to the exponents are not expected to erase the large 
gap z ~ 2/3. 
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4. 2-POINT FUSION RULE FOR 4-POINT 
CORRELATION FUNCTIONS: THE EQUATION 
FOR THE EXPONENT p 


In this section we use the fact that we know exactly 
the ultraviolet exponent A to derive the fusion rule for 4- 
point correlations. One immediate consequen ce o f these 
rules will be a derivation of the scaling law (1.9). The 
n-point simultaneous correlation function of Eulerian ve¬ 
locity differences is defined as 


F“^-“(R 1 ,R 2 ...R n ) 

= ( 6 u a (r a \R. 1 ,t) 5 up(ro\R. 2 ,t)... 5u u (r 0 \Tl n A)) . (4.1) 

where (5u was defined by ( |2.2| ) such that Rj = Tj — r 0 . 
By “fusion rules” we mean the scaling structure of F n 
when two or more of the vector separations Rj tend to 
zero or to each other. In this section we discuss 2-point 
fusion rules. Without loss of generality we can denote the 
coalescing coordinates as ri and r 2 , and we will consider 
all distances of all the other coordinates r 3 ... r n from ro 
to be of the same order 0(R) 7'i,r 2 . To understand 

the fusion rules we will use diagrammatic language; at 
the end of the discussion one can recast the results into 
formal operator algebra of a new type without reference 
to diagrams. This operator algebra (which allows nrulti- 
scaling) is an exciting subject that will be taken up fully 
in a separate publication |lGj . 


A. Classification of the Diagrams for the 4-Point 
Correlator 


diagrammatic presentation of this object is shown in Fig. 
|l] panels c and d. The right part of diagram 1 is the three 
point object A shown in panel b of Fig. in terms of the 
dressed vertices A,B and C, and the dressed propagators. 
The left part of diagram 2 in Fig. is the third order 
correlation function F 3 . Its diagrammatic representation 
in terms of dressed vertices and propagators is given bv 
Fig.7 of paper II, and in a more compact form in Fig. |Gj 
panel c. We find again in this case the same three-point 
objects A and V that appeared already in panel b of this 
figure and in panel d of Fig. |T]. 

In addition to these weakly linked diagrams we have 
the strongly linked diagrams appearing in Fig. |j] panel b. 
All these are diagrams whose resummed parts are con¬ 
nected via two-propagator bridge. The strongly linked 
diagrams are constructed as follows: locate the princi¬ 
pal cross section of a diagram in the infinite expansion 
of F 4 . Move from the principal cross section to the right 
and to the left until you find the first two-eddy reducible 
link, which is going to form the two-propagator bridge. 
The object between these two bridges is a contribution 
towards one of the resummed central objects on the RHS 
of the equation in Fig. [?], and see Fig. || for the begin¬ 
ning of the diagrammatic expansion for these objects. At 
this point collect all the diagrams with the same central 
object and resum them to the left and to the right. The 
result of this resummation is shown in Fig. |?]. The four 
point objects on the right and on the left of the central 
objects are exactly the same as those appearing in Fig. 

W 


B. Asymptotics 


In this section we will analyze the 4-point correlator 
and the dissipation correlation fun ctio n K ee . In particu¬ 
lar we derive the scaling relation (0 for the exponent 
that governs the power law decay of K ee according to 
the (ILq). The derivation begins with the examination of 
the 4-point correlator ... R 4 ). This quantity is 

represented as an object with 4 wavy tails, each of which 
represents one of the velocity differences with coordinates 
r 4 ... r 4 . The trivial contributions to F 4 , which we de¬ 
note as F ^ 0 " 1 , are the three different products of F 2 which 
graphically are the “unlinked” (or “reducible”) diagrams 
which stem from the Gaussian decomposition. The next 
set of diagrams that contributes to the 4-point correlator 
is shown in Fig. |j. These are all the weakly linked dia¬ 
grams in which the coordinates X\,X 2 are linked to the 
coordinates x 3 , X 4 via one propagator. As in the case of 
Gn,m we cannot have repeated one-propagator bridges 
since such contributions are resummed in the Dyson- 
Wyld line resummation. There are two possible types 
of bridges: (i) the bridge ends with a straight line and 
(ii) the bridge ends with a wavy line. In case (i) the di¬ 
agrams can be resummed into diagram 1 of Fig. |] panel 
a. In case (ii) the diagrams resum into diagram 2. The 
left part of diagram 1 is the Green’s function G 2.1 ■ The 


1. Strongly Linked Diagrams 

Assume that r 4 , r 2 ~ r are much smaller than 
7 * 3 , 7*4 ~ R but still in the inertial interval. Denote 
the structure whose outer coordinates are r a , rf,, r 3 , r 4 as 
^m(r a ,r6, r 3 , r 4 ) with m = 1... 4 according to the enu¬ 
meration of the diagrams in Fig. [?J. Diagram (1-4) can 
now be written analytically as 


(1) -J 

[ dXadx b g 2 ,2(jco\x 1 ,X2,Xa,x b )'it 1 (r a ,r b ,r3,r4) , 


(4.2) 

(2) u 

[ dx a dx b g 3jl (r Q \x 1 ,X2,Xa,x b )'i>2(ra,'rb,r3A4) , 


(4.3) 

M-J 

[ dx a dx b Q 2 , 2 ( r ol^i, £ 2 , £aj Xf,)H/ 3 (r a , rb, r 3 , r 4 ) , 


(4.4) 

(4) u 

[ dx a dx b g 3j i(ro\xi, x 2 , x a , x b )^4(r a ,Tb, r 3 , r 4 ) . 


(4.5) 
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Suppose that we can argue that all the functions 
'I' m (r a , rh, r 3 , r 4 ) contribute at r a , V b sm aller contribu¬ 
tions than those of <h(r a ,rf,) in Eq. ( 3.15| ). If so, we can 
think of these functions as generalized perturbations 
SD that leave the n 2 dependence universal. According 
to ( 3.12| ) and ( 3.13| ) the main n 2 dependence comes from 
diagrams (1) and (3) and it is r^. 

To demonstrate that this is indeed the case we can look 
at a typical diagrammatic contribution to any of the di¬ 
agrams in Fig. [i]. For example, the first contribution 
to 4'i(r a ,!•{,,r 3 ,r 4 ) is shown in Fig. |[ This particular 
example comes from diagram (1) in Fig.q, panel b upon 
taking diagram (1) for the rung in Fig. |. Take r Q ~ r\, 
and examine the dependence on r a . The two vertices 
contribute r” 2 . The correlator in between is worth r’ 2 . 
The two stretched correlators connecting to larger coor¬ 
dinates contribute r 2< ’ 2 . The two time integrals in the 
vertices c and d contribute r 2z where z is the dynamic 
scaling exponent, r(r) oc r z (and see details in paper I). 
In total we have rf with (3 = 3(^2 + 2z — 2. This has to 
be compared with the expo nent a of <I>(r a ,rf,) which is 
a = (4 —2 according to &■ In K41 evaluation [3 = 4/3 
whereas a = —2/3. The gap is so large that small correc¬ 
tions to K41 cannot change the fact that the exponent of 
$ is dominant, as we want. Note that this demonstration 
is perturbative (order by order with only the lowest order 
done explicitly), but again the magnitude of the gap is 
sufficient to validate the claim. 

Ex actly the same argum ents hold f or th e contribution 
(4.4). The contributions ( |4.3| ) and (15) relate to the 
asymptotic behaviour of the other Green’s fun ction G 3 .i- 
It is clear however that due to the result in (3.2C) the 
exponents are the same. 


2. Weakly Linked, Diagrams 

The analysis of the weakly linked diagrams in this case 
follows closely the discussion of section 3C. The diagrams 
can be resummed as shown in Fig|], panel b. We have 
two contributions; in diagram (a) we have on the left F 3 
and on the right the object defined in the second line of 
the figure. In diagram (2) we have Q 2 ,i (shown in Fig. 
0 ) on the left and the object defined in the third line of 
Fig. |], panel b. In the limit r±, r 2 ~ r <C r$, r 4 ~ R the 
integrals over rt> in both diagrams contributes mainly in 
the region r t ~ R. Accordingly in diagram (2) we have 
the same situation that was discussed in section 3C, and 
see Fig||, panel b. Diagram (1) on the other hand may be 
bounded from above by the asymptotics of F 3 with one 
large and two small coordinates. Under the property of 
rigidity one can see that this correlator is independent of 
the large coordinate, and is therefore proportional to r 1 . 
Again we have a gap, and this gap can only increase if we 
take into account the time integral of the Xb vertex. Thus 
the weakly linked diagrams can be again disregarded. 


C. 2-Point Correlation of the Energy Dissipation 


At this point we ca n em ploy our results to evaluate the 
correlation function ( |l.4| ) . We first write this quantity in 
a way that makes its relation to F 4 clear: 

K ee (R) = v 2 lim lim Vi ct V 2 / 3 V 3 7 V 45 

ri2,r3 4 —>0 ri 3 —>R 

pa ft 1 5 ( ri)I . 2)I . 3;r4 ) _ F a,/3 '(ri,r2)F 7 ' 5 '(r3,r4) 


$a[3$a'0' T ^ad'^a' 


aft' Oa' (3 


T Sryb' dry! b 


(4.6) 


This is a generalization of the situation discussed above 
in the sense that we have two pairs of coalescing points, 
ri,r 2 and r 3 ,r 4 . By examining Fig. 8 a we can see that 
the symmetry allows us to treat each coalescing pair sepa¬ 
rately. In the limit ri 2 and r 34 being much smaller than R 
each pair of derivatives ViV 2 and V3V4 will contribute 
a divergent term proportional to rj~ 2 A and r^" 4 A respec¬ 
tively. The dependence on R can be found knowing that 
the overall exponent of F 4 is 1 / 4 . Thus the over all scaling 
exponent of the quantity in the RHS (4.6) must be (4 — 4. 
Using the fact that A = 2 — £ 2 we conclude that 


ViV 2 V 3 V 4 F4(n,r 2 , r 3 , r 4 ) 

i [ri 2 r 3 4] _A i? C4 “ 2C2 ^ 4/3_C4 


(4.7) 


e*/3r 


where for simplicity we suppressed the vector indices and 
we made the relation dimensionally correct by introduc¬ 
ing some renormalization scale £. Note that for K41 scal¬ 
ing exponents the factor e 4//3 fixes the correct dimension¬ 
ality of the RHS. For anomalous exponents one needs 
a renormalization scale to take care of the difference be¬ 
tween the K41 and the actual value of the exponent £4. It 
will be shown later that this renormalization scale must 
the outer scale L. At this point it do es n ot matter. The 
divergence in the limit indicated in (4.6) should be un¬ 
derstood in light of the full theory for F 4 , in which the 
I'-difFusive terms are explicit. The role of these terms is 
prec isely to truncate the divergence that is implied by 
(4.6). As a consequence the divergence is only applicable 
in the inertial range with ri 2 , r 3 4 > 77 , whereas in the 
dissipative regime the divergence disappears. Thus for 
evaluating K ee (R) via inertial range values we must re¬ 
place the limit ri 2 ,r 34 —> 0 by ri 2 = r 34 = 77 . Thus we 
can write 


K ee {R) ~ g4/3 r? -2A i? C4-2C 2 ^4/3-C4 _ ( 4 . 8 ) 

Finally, we need to evaluate the viscous scale r /. This is 
done from the definition of e which is 

e ~ v lim ViV 2 S 2 (ri 2 ) . (4.9) 

ri2—>ri 

In the same spirit we can write the structure function 
5 2 (i?) as 


S 2 (R) ~ e 2/3 R C2 l 2/3 -« 2 , (4.10) 
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where again the renormalization scale l fixes the dimen¬ 
sions. From the last two equations one can compute 


?V3 


-A£2/3-C2 _ 


Substituting this in (L 8 ) we get finally 

2C2—Ci 


_2 [ i\ 


Kee(R) ~ r [ -J 


(4.11) 


(4.12) 


which is the scaling law (L 8 ) that we wanted to derive. 
It will turn out that the renormalization scale i i s th e 
outer scale of turbulence L. From this we get Eq. (Ill) 


H = 2C2 - C4 • 


(4.13) 


5. 2-POINT FUSION RULE FOR MANY-POINT 
CORRELATION FUNCTIONS: THE 
DERIVATION OF J N 

In this section we evalu ate the scaling exponent of the 
dissipative terms J n ( ^2.40 ) of the balance equations. The 
strategy is to expose the divergence with the ultraviolet 
cutoff 77 for which we have an exact evaluation of the 
exponent A. Together with the overall scaling which is 
determined by ( n we will be able to compute the R de¬ 
pendence of J n (R). Having that we will use an argument 
due to Kraichnan to find the coefficient of the power law. 
We believe that our result for J n (R) is exact, even though 
some of steps taken are not rigorous. Afterwards we will 
evaluate the interaction term D n ( R ) and use the balance 
equation as a non-perturbative constraint to deduce the 
scaling exponents 


A. Classification of the Diagrams as Weakly Linked 
or Strongly Linked 

The n-point correlation function F n can be represented 
symbolically as an object with n wavy tails, each one 
representing one of the n coordinates r j associated with 
a velocity difference i5u(ro|Rj). We remind the reader 
that for simultaneous correlation functions one can use 


either the Eulerian differences (T2) or the BL-velocity 
differences ( 2.1 1| ) , since they are the same at t = 0. In the 
diagrammatic expansion for the simultaneous F n there 
appear time dependent correlations, and there the theory 
calls for the use of BL-velocity differences. 

Having two special coordinates r a and we can ask 
how the part of the diagram containing these coordinates 
is linked to the rest of the diagram. This part can be 
connected via one propagator, see Fig. E panel a, via 
two propagators, see Fig. [it] panel b, or via three or 
more. As in the case of the 4-point correlator and 4-point 
Green’s function a one-propagator bridge cannot appear 
again between the legs carrying the designation a, b and 


the body of the diagram. In total one has n(n — l)/2 
weakly-linked contributions in each of which the role of 
the weakly linked pair is played by one of the available 
pairs of legs. In addition one has also double weakly- 
linked contributions with two bridges made of a single 
propagator which connects different pairs, etc, and see 
diagram 3 in Fig. |l(] panel a. Such diagrams do not 
play an important role in our analysis. In displaying the 
diagrams in Fig. 10a we have included the bridge (which 
is a correlator or a Green’s function) in the object on 
the left of the bridge. For this reason the object on the 
right of the bridge begins with a vertex x c . According to 
our convention this vertex is denoted by a small empty 
circle. The diagrams appearing in Fig. [H] panel a should 
already look familiar. In diagram 1 we again have a F 3 
correlator on the left, integrated over x c with an (n— 1 )- 
point object on the right. This is the generalization of 
the diagram 2 in Fig. [(] panel a. In dia gram 2 we have 
the 3-point Green’s function C/ 2,1 of Eq.(3T) on the left, 
again integrated over x c with a different (n — l)-point 
object. This can be compared with diagram 1 in Fig. |(i| 
panel a. For our considerations the precise nature of the 
object on the right is irrelevant. 

Similarly, all the contributions with two propagators 
serving as links can be resummed into the objects shown 
in Fig. [loj panel b. Again we included the two prop¬ 
agators that form the bridges in the objects appearing 
to the left of the bridges. Again this results in the n- 
point objects on the right being attached to the left by 
two vertices x a and X/,. Their interpretation is as fol¬ 
lows: diagram 1 is a fourth order correlator linked via 
integrals over x a and Xb to the rest of the diagram. In 
diagram 2 the linked object is Qz,i, (3.4). In diagram 3 
the linked object is precisely our familiar second order 
Green’s function C/ 2 , 2 • 

Links with three or more propagators have been al¬ 
ready taken into account in this presentation. This classi¬ 
fication of the diagrams is based on starting with the two 
special coordinates X\,X 2 on the left, and then moving 
to the right and stopping when the first one-propagator 
bridge appears. All these diagrams are resummed ex¬ 
actly into one of the contributions in Fig.10a. If there 
is no one propagator bridge, we start again from the left 
and monitor all the two-propagator bridges, ending with 
the last pair. All such diagrams are resummed into one 
of the contributions in Fig. fill, panel b. 


B. Asymptotic Tensor Structure and Fusion Rules 

In this section we find the tensor structure of F^'" u 
when the two coordinates r\, r 2 are much smaller than 
the rest. When all the coordinates r, 3 ,... r n are of the 
order of R r, whereas 7~i, 7*2 are small, the property 
of rigidity that was demonstrated in paper II requires 
that the main contribution to the integrations over r a 
and Tb in the diagrams of Fig. 10 panel b come from the 
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region r a ~ r b ~ i?. In this case we have a very similar 
situation to the one discussed before in the context of 
Fig.0 panel b. Accordingly the three objects on the left 
of the double bridge of the diagrams in Fig. [l(] panel b 
are familiar and have the same scaling with respect to r, 
i.e. r^ 2 . The weakly linked contributions shown in Fig. 
|TTi| panel a have the same objects on the left of the bridge 
as those shown in Fig. |(i| panel a for F&, and they are 
irrelevant for the same reasons. 

Limiting our attention to strongly-linked diagrams 
with two propagator links we examine now the tensor 
s truct ure of We need to keep in mind equation 

( |3.13 ) which means that when r\,V 2 —» 0 5F a P oc F a P. 
In light of Fig. || the objects and f 7<5 have the 

same asymptotics. We have found before that they are 
proportional to F a!3 {r\,V 2 ). Diagram 1 is again propor¬ 
tional to the same quantity in light of Fig. [?] panel b. We 
can therefore write the following fusion rule: 


lim F^ 7 -“(r 1 ,r 2 ,R 3 ...R 7l ) 

rL,r*2—>0 

= F^(( r 1 ,r 2 )^(R 3 ...R n ) • (5.1) 


general forms for isotropic conditions. Incompressibility 
dictates the value of a m but we do not need to compute 
it for our purposes. Lastly we note that 

F ary (dpep, dpep) = S^{dp) oc df 

(5.6) 

Incompressibility requires the relation ai = ( 2/(1 + ( 2 )- 
Pre sently we can substitute all this knowledge into Eq. 
m - One should note that in the inertial range the ve¬ 
locity field is not smooth, (£2 < 2 ) and we may run into 
the danger that the quantity computed depends on the 
angle between the vectors d, R. However the procedure 
implied requires taking the limit such that all dp = r). 
This can be checked explicitly by introducing the ten¬ 
sor structure to all the quantities appearing in the limit 
that indeed the condition that all the components dp are 
the same guarantees that the limit is independent of the 
angle. The result of the substitution is 

J2m{R) = 2 mC2m^-^R^-^n 2rn/ ^-^ . (5.7) 


Rq,R- y 

Ocx.'y 


R 2 


Here ’P „_2 is a homogeneous function of its arguments 
when they are all in the inertial range. The scaling expo¬ 
nent of \kn _2 is C, n — C.2- The reason for that is clear: the 
scaling exponent of F n is £ n , but one £2 is already carried 
by F a/3 . The tensor structure of T 7 !'^ is not known in 
detail at this point, except that it conforms with incom¬ 
pressibility and isotropy. In other words, this quantity is 
independent of the vectors ri,r 2 - 


C. The Evaluation of J„(R) 


At this point we are ready to evaluate the dissipative 
ter ms in the balance equation. The equation to consider 
is ( 2.40| ). As was done in the context of the evaluation of 
K ee we evaluate the quantity in the inertial interval, and 
then replace the limits dp —> 0 with the ultraviolet cutoff 
dp = 77 . In the inertial interval we write 


T£%(dpep,dpe b eta,R) = F aa {dpep,dpep)'i! < ^_ 2 {R) , 

(5.2) 

T^idpep, dpep, R) = F“ 7 (d / 3 e /3 , d^e^^R), 

(5.3) 

where the fusion rule (H has been used. The two func¬ 
tions 4/ satisfy 


*&- 2 (fl) = ^ am ~ Ca 


- 2m—2 

*^(R) = 


1-C2 


R a Rj 

Oryy CLq 


R 2 


(5.4) 

(5.5) 


with A\ and A 2 being constants. The reason for these 

forms stems again from the fusion rule. T- 1 ) and 

are scalar and 2 -tensor respectively, and we wrote their 


We again use a renormalization scale t to fix the dimen¬ 
sions. It will turn out that this renormalization scale is 
the outer scale of turbulence L. 

We will see that this form of J n is sufficient for the 
calculation of the exponents C, n through the use of the 
balance equation only if the dependence of ( n on n is 
linear (i.e in the /3-model). If the dependence is nonlinear 
(multiscaling) we need to determine t he c oefficients C 2m 
exactly. To do so we rewrite now Eq. ([.5 . 7| ) in terms of the 
structure functions S n . This way of writing is compelling 
only when the scaling exponents are a nonlinear function 
of n, as will be clear in a moment. We write S n in the 
form 


S 2m {R) ~ e 3m / 3 R^™f m / 3 -^ 


(5.8) 


This form is the generalization of ( 4.10| ) and is the most 
general form that conforms with scaling and is dimen¬ 
sionally correct. Using Eqs. (4.1C) and ( [4.11 ) in (5.8) we 
find the convenient representation 


T C 2m S 2m (R) 

J2m{R) = m-^—J 2 - 


C 2 


S 2 (R ) 


(5.9) 


We stress that this result is valid only when R 77 , since 
we used the asymptotics, and only when is nonlinear 
in n. If the scaling exponents are linear in n we can 
have other contributions like S 2m +i/ S3 or any othe r ra tio 
whose scaling exponent is C, 2m — C, 2 . In particular (T9) is 
not applicable to Burgers turbulence. 

Finally, we will employ an idea that is due to R.H. 
Kraichnan |l^| to argue that in the multiscaling ca se the 
coefficient C 2m is ?n-independent. Begin with Eq. ( 2.22| ) 
which is rewritten as an integral over the distribution 
function P(w): 
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J2m(R) = j dv/P(yv)w 2{ - m 1} l«a([V? + V?> a |w) . 

(5.10) 

Here ([Vj! + Vji ]ty a |w) is the conditional average of 
[Vj: + V 2 0 ]w a conditioned on a given value of w(r 0 |r 0 + 
R, t ). The point to observe now is that the only way to 
recover our result ( |5.9| ) when £ n is a nonlinear function 
of n is to demand that the conditional average satisfies 


/ r 2 2 , . \ w Q (r 0 |r 0 + R,t) 

\[V r + V r „]u> a |wy = C -- , (5.11) 


where C is some coefficient which is evidently indepen¬ 
dent of m. It follows that C 2m is independent of to. 

Note that this result for the conditional average is only 
valid in the inertial range of scales, since it has been de¬ 
rived using Eq. (^1]) whic h is o nly valid there. Notwith¬ 
standing we can employ ( |5.11| ) right aw ay to compute 
the vector quantity J 2 rn +i(R)- Writing ( |2.27| ) again as 
an integr al ove r the distribution function R(w) and sub¬ 
stituting (5.11) we find 


J. 


2 m+l 


(R) = 


(2 in + 1) ^2m+l(R) 

2 2 S 2 (R) 


(5.12) 


We can thus summarize this section with a result that is 
valid for both odd and even n by using the scalar coun¬ 
terpart of the vector quantities: 


Jn{R) — J2 


nS n (R) 

2S 2 {R) 


(5.13) 


This is the final result of this section. We note that such 
a scaling formula for J 2m was suggested by Kraichnan in 
the context of passive scalar advection ( 0 ) , and was 
derived in |Ig| |. 


D. The Dynamical Exponent z„ 


z n = (,2 for all n . (5.15) 

This is a non-trivial prediction that to our knowledge has 
never been tested either in experiments or simulations. It 
appears to be exact. It is interesting to notice that inde¬ 
pendently of the question of multiscaling in the spatial 
scale, the temporal scaling is simple. 

It is amusing to try to understand (|5.15| ) intuitively. In 
doing so we want to separately understand why z 2 = C, 2 
and then why all z n are the same. The first finding 
seems to contradict the naive dimensional evaluation of 
t 2 (R) as R/^/S 2 (R), which is the “turn-over” time of 
R-eddies with characteristic velocity \JS 2 (R). This eval¬ 
uation would lead to z 2 = 1 — ( 2/2 which is wrong. 

Another way of thinking that leads to the right result 
is to estimate the rate of energy dissipation as the ra¬ 
tio of energy of R-motions, which is S 2 (R ), by the time 
scale r 2 (R). Since the rate of energy dissipation is R- 
independent (being e), this fixes t 2 (R) to scale as R^ 2 . 

The n-independence of z n is more subtle, and we posto- 
pone its discussion to the forthcoming paper [pT|. Here 
we just want to point out that this result entails a pre¬ 
diction about the measurement of dimensionless ratios of 
structure functions, like S^/S^ 3 , Si/Sf etc- ' n decay¬ 
ing turbulence. The prediction is that such relatios are 
R-dependent but not time dependent. We believe that 
this is not in contradiction with what is known about 
decaying turbulence behind a grid. 


6. THE INTERACTION TERM IN THE 
BALANCE EQUATION 


In this section we pres ent the anal ysis of the interaction 
term D„, and see ( 2.3(i| ) and ( 2.37 ). The question that 
was left at the end of sec 2C 1 is whether the integral 
over ri converges. Order by order analysis of the type 
presented in paper I indicates that the answer is yes. 
However we need now to consider the nonperturbative 
answer using what we have learned so far. 


As a windfall profit from the calculation of J n we can 
show now that the dynamical exponents associated with 
the n-point correlations are all the same. The dynamical 
exponent z n is defined by the assumption of scale invari¬ 
ance of the n-point, n-time correlation function of the 
BL-velocity differences in the form 

F n { ARi, A ... XR n ,\ Zn t n ) 

= A c "E n (Ri,fi,...R„,£ n ) . (5-14) 

The meaning of this form is that the typical time 
scale associated with n-point quantities scales like R Zn 
when R is changed. In particular we can exa mine E q. 
( |2.20| ) and realize that J n , which according to (B is 
nJ 2 S n (R)/S 2 (R), gives us the desired time scale. In fact, 
the R dependence of the time scale is determined entirely 
by S 2 (R), and therefore 


A. Locality of the integral in the interaction term 

In order to do this we need further asymptotic prop¬ 
erties of the functions T n which appear in the inte¬ 
gral. For brevity we will suppress the tensor indices of 
these objects, and to consider even and odd n in the 
same way. The convergence of the integrals depend on 
T n (R+ri, R+ri, R) and T n (ri, ri, R) when n R and 
when ri R. So far we have only analyzed T„(ri, ri, R) 
when ri <C R. The full analysis of the two unknown 
asymptotics is as involved as the one presented above, 
and we will present them in a separate publication [fni|| . 
Here we will simply employ the results that we need for 
the present analysis. 

Consider first T n (R + ri,R + ri,R) for r\ small. The 
analysis in jn| shows that for ri <R 
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T„(R + ri,R + ri,R) - S n (R) ex S 2 {ri) oc r [ 2 . (6.1) 

Next consider T„(ri,ri,R) in the limit r± R. The 
analysis in [HI leads to 


T n (ri,ri,R) ex R c 


1-2^ Cn-2 

' 1 


( 6 . 2 ) 


These results can be employed now in the integral for 
D 2m . In this integral we have the projection operator 

P, which has a delta function and a longitudinal part. 
It was demonstrated in section 2 that the delta func¬ 
tion leads to the expression ( [2.34 ), as if there were no 

pressure. The longitudinal part of P(ri) is proportional 

to 1/rf. The integral J driP(ri) by itself diverges log¬ 
arithmically. The rest of the integrand (i.e. dT n /dr i 7 ) 
behaves like rf A r i 7 . Simple power counting indicates 
that the integral diverges on the whole in the ultraviolet 
region. In fact, this power counting is misleading, since 
the integration over the angles vanishes. The projection 
operator is an even function under the inversion of ri, 
whereas the leading term of the rest of the integrand is 
odd. The next term in the expansion of dT n /dr i 7 is even 
under the inversion of ri, and is of the order of rf 2 . The 
resulting integral f drirf 2-1 converges in the ultraviolet. 

Note that this analysis indicates that each of the two 
terms in the integral for D n converge in the ultraviolet 
independently. In fact, we see from Eqs. dj and ( |5.2| ) 
that the two terms have precisely the same asymptotics, 
and they ma y exa ctly cancel in the limit. In addition 
we see from ( 2.36 ) that the leading asymptotics of the 
two terms cancels exactly also in the infrared. These 
facts are important; we will argue below that the leading 
scaling behaviour for D n which is naively calculated from 
each term separately may cancel, and the actual scaling 
is determined by the next order contribution. This will 
be a mechanism for multiscaling. 

Notwithstanding the exact cancellation of the lead¬ 
ing infrared behavior we need to examine the infrared 
convergence of the integral. Each one of the terms 
in the integrand of D 2m has the asymptotic form 
J(dri/ri)dr < [ r ‘ ( ’ 7 l ~ 2 /dr i which converges separately in 
the infrared. The difference should converge even faster. 

In summary, we argued here that the proof of local¬ 
ity of the integral for D 2m extends beyond order by or¬ 
der considerations. Similar arguments allow reaching the 
same conclusion for the integrals in Z^m+i- 


B. Rough Evaluation of D n 

The conclusion of the last subsection is that the main 
contribution to the integrals appearing in D n comes from 
the region rq ~ R. Accordingly the integral can be eval¬ 
uated as discussed in [fl5l: 


D n (R ) 


dS n+ 1 (R) 

dR 


(6.3) 


This is exactly of the form computed in Eq. (|]||) for 
the case without pressure. We will see next that if this 
evaluation is to be trusted, its unavoidable consequence 
is that the scaling exponents C, n are a linear function of n. 
We will see that this allows only K41 scaling o r /3- model 
type scaling It is possible however that ( |6.3| ) is an 
overestimate; we ar gued above that the interaction term 
as shown in Eqs. ( 2.36 ) and ( 2.37 ) may have a cancel¬ 
lation of the leading scaling behavior which is valid for 
every one of the terms in the integrals separately. We will 
therefore also study now the next order term that will be 
the proper evaluation of D n if the leading evaluation in¬ 
deed cancels. We will see that the resulting evaluation 
culminates in multiscaling in close agreement with ex¬ 
perimental observations. Moreover, we will argue that as 
far as we can see within the new approach developed in 
this series of papers, the scenario that involves the can¬ 
cellation of the leading scaling behaviour of D n is the 
only scenario that allows multiscaling in turbulence. It is 
possible that there exists a symmetry or a sum rule that 
leads to such a cancellation but we do not have at the 
present time a theoretical proof or even a good argument 
as to why and how it happens. 

To evaluate the next order scaling contribution of D n 
we need to return to the diagrammatic expansion of F n , 
Fig. [h]. In the discussion in section 5 we explained 
that in the asymptotic regime of two small coordinates 
the weakly linked diagrams are negligible compared with 
the two-propagator bridged contributions, even in K41 
scaling. Now however we are interested in these dia¬ 
grams when all the coordinates are of the same order, 
and it is evident that K41-wise they all have the same 
scaling with R. In fact, the unlinked contributions to 
F n which are obtained from the Gaussian decomposition 
(i.e. all the contributions F p F q with p+ q = n) also have 
the same K41 evaluation. In addition we have a set of 
weakly linked contributions such as the ones displayed in 
Fig. nq, panel a. Again they have the same scaling in 
K41. Accordingly we need to think which contributions 
are dominant when the leading scaling (. R/Lp n cancels 
in the evaluation of D n . (Since we are going to show 
that anomalous scaling requires the normalization scale 
to be A, we assume this in the present discussion without 
further ado). 

The estimate of the scaling exponents of all these var¬ 
ious contributions is facilitated by the fact that we are 
interested now in the “local” situation when all the co¬ 
ordinates are of the order of R. Thus for example the 
diagrams 2 in Fig. ^ panel a are redrawn in Fig. |nj. 
The objects in the left grey ellipse in diagrams 1 and 2 
are representation of F 3 shown in panel c of |(]. The ob¬ 
ject in the right ellipse of diagram 2 belongs to F n _ 1 . 
One can see this by taking n = 3 and looking back to 
panel c of Fig. £s| However in diagram 2 we counted the 
F '2 bridge twice. The overall scaling exponent is there¬ 
fore £3 + £ n _ 1 — ( 2 - Using Fig. |(] panel b it can be also 
seen that the diagram 1 in Fig. |ll] has the same scal¬ 
ing exponent. Similarly we can analyze the diagram 1 in 
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(7.2) 


Fig. [lC] panel a and argue that its scaling exponent is 
C 3 + Cn—i - <2. etc. 

Now we need to understand which of these contribu¬ 
tions will take the lead if the main scaling (R/L)^ n can¬ 
cels. To guide our thinking we will assume that the scal¬ 
ing exponents are neither K41 nor /3-model, but are non¬ 
linear functions of n. Holder inequalities then require 
that the increments between and C™-i will be nonin¬ 
creasing functions of n , i.e. 


S, 


Cn+l C n — Cn Cn —1 


(6.4) 


With these inequalities one sees that the unlinked contri¬ 
butions which scale with (i?/i)G+G are always smaller 
than the weakly linked contributions, and that of all the 
weakly linked contributions the leading one is the one 
which we singled out in Fig. ^ panel a, with the scaling 
exponent £3 + (n-i - ( 2 - 

The meaning of this result is that instead of evaluating 
T n in the integrals for D n as S n (R ) we need to evaluate 
it as S n -i(R)S 3 (R)/S 2 (R). Correspondingly the evalua¬ 
tion (|6.3|) changes to 


D n (R) = d n { Cn) 


S n (R.)S 3 (R ) 
RS 2 (R) 


(6.5) 


where d n ( Cn) is a coefficent which depends on the numer¬ 
ical value of the scaling exponent. 


7. ANALYSIS OF THE BALANCE EQUATION AS 
A NONPERTURBATIVE CONSTRAINT 

At this happy moment we can use all the knowledge 
accum ulate d so far to go back to the balance equations 
([2.24 ) and (2.25) that we rewrite in the form 


D n (R) = J n (R) 


(7.1) 


where the forcing term is not displayed because it is neg- 
ligi ble, c f. section 2C3. The evaluation of J„(_R) is give n 
in (5.13). The evaluation of D n is either ([k^) or ( |6.5| ), 
depending whether there is oris not cancellation in the 
leading scaling behaviur of D n . We will show now that 
option (|6.3| ) leads inevitably to linear scaling, whereas 
option (p.5|) leads to multiscaling. 


A. Scenario of Linear Scaling: Burgers Turbulence 
and the /3-model 


The evaluation (|6.3D is exactly correct only when there 
is no pressure term and the projection operator is a delta 
function. This is the situation for example in Burgers 
turbulence f 2 C|| . It may or may not be a proper evaluation 
of D n also in the case of Navier-Stokes turbulence, as 
discussed above. We examine now the consequences of 
this evaluatio n wh en su bstitu ted in t he balance equation. 
Substituting (6.3) and ( 5.13 ) in ( ]7. lj ) we find 


n +i(R) _ .S n (R) 


R 


s 2 {r ) ■ 


For n = 2 we recover the known result that S 3 (R) ~ eR. 
Accordingly we can rewrite (7.2) as 


S n+1 {R)S 2 (R) ~ S n (R)S 3 (R) . (7.3) 

In terms of the scaling exponents this results reads 


Cri+l + C2 ~ Cn + C 3 • 


(7.4) 


The only solution of this equation is the linear law Cn = 
a + bn , where a and b are some consta nts. Knowing 
that ( 3 = 1 and using our scaling law ( 0 ) which is 
H = 2( 2 — £4 we find that the only solution is 


A n (n — 3) 

U - g - n 3 


(7.5) 


It is interesting to note that this result is identical to the 
prediction of the /3-model |li|, coefficients and all. This 
should not surprise us too much. After all, once we have a 
linear dependence two constraints fix the linear law com¬ 
pletely. Note that this law includes as a special case the 
exponents of Burgers turbulen ce w hich are Cn = 1 for all 
n (2^]. This is obtained from (7.5) when ft = 1. Never¬ 
theless the full analysis of the Burgers equation using the 
techniques developed in this series of papers needs special 
attention due to the importance of the incompressibility 
constraint in so many of our calculations. The Burgers 
case deviates so strongly from K41 that the issues of lo¬ 
cality and rigidity of the various diagrams needs to be 
assessed separately [^l| . 

We can also show now that CD implies that the 
renormalization scale is th e ou ter scale of turbulence L 
as claimed before. Eq. CD was derived by assert¬ 
ing that the Gaussian contribution to J n is negligible 
compar ed to the connected ladder contributions which 
led to ( 5.13| ). The Gaussian contributions are domi¬ 
nated by eS n -2 whose scaling exponent is £ 2 . If we used 
this c ontr ibution as the leading one in the balance equa¬ 
tion (0) we would have obtained the scaling relation 
Cn+i = (3 + Cn -2 with the obvious boundary condition 
£0 = 0. The solution of this recursion relation is K41 
scaling with ( n = n/ 3. If this is to be rejected in favor of 
(7.5) the Gaussian contributions must be indeed smaller 
than the ones we kept. The conclusion is that 


S 2 {R)S n - 2 (R) < S n (R) . 

In turn this inequality implies that 


R 


Cn -2 


< 


R 


Cn -<2 


(7.6) 


(7.7) 


Substituting the scaling exponents from (|7.5|) we con¬ 
clude that ( R/iY < 1, which can only happen if R < t 
for any R in the inertial inteval. This identifies t with L. 
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B. Scenario of Multiscaling 


Substituting 


and (5.13) in © we find 


S n (R)S 3 (R) nJ 2 S n (R) _ nS n {R)S 3 (R) 
n[U> RS 2 (R) 2 S 2 (R) 2 RS 2 (R) ’ 

(7.8) 


where the last form is obtained from J 2 = S 3 /R. From 
the point of view of scaling exponents this equation is an 
identity. The only way to compute the exponents </ n now 
is from the coefficients in the balance equation: 


2d n (()n) — tl • 


(7.9) 


To this aim the rough evaluation of D n in section 6 B is 
not sufficient; we need to be much more precise in order 
to compute the („ dependence of d n . 

Clearly, the computation of coefficients is exceedingly 
hard. Only when we have exact form for the functional 
dependence of the many point functions we can hope to 
compute the coefficient. We did have an exact form for 
J n because we understood how to resum its diagrams; 
consequently we believe that we have computed the coef¬ 
ficient of J n exactly. D n is a different matter; at present 
we do not have an exact functional form for it. Order by 
order considerations are not helpful for this issue, and we 
still do not know how to exactly resum the diagrammatics 
for D n . We will therefore try to guess the £„ dependence 
of the coefficient of D n . 


hiR R2 


S 3 {R)R d 
S 2 {R) dR 


S n (R) , 


(7.11) 


with b being now an n-independent coefficient. The phys¬ 
ical meaning of this guess is that the R-dependent eddy 
viscosity h(R) takes here the form 

HR)=b Wr- (7 - 12) 

Note that the eddy diffusivity which is introduced here 
scales like R 2 ~= R A . 


2. Scaling Exponents in the Eddy-Viscosity Approximation 


Using ([7.11) we compute 


<*»««) = &Cn(3 + Cn-C2) , (7.13) 

where we used the fact that £3 = 1 . Together with © 
we find a quadratic equation for ( n : 


2b(n(3 + (n - ( 2 ) =n . (7.14) 


We remind the reader that this equation is valid for n > 2 
since there is no cancellation of the leading scaling order 
in D 2 . Using again the fact that <^3 = 1 we find the 
constraint 


26(4 - C 2 ) = 3 


(7.15) 


This leaves us in (7.14) with only one unknown number 
which we take as Q 2 . Solving for Q n we find 


1. The Eddy-Viscosity Approximation of D„(R) 


To guide our thinking we recall some results from the 
theory of passive scalar advection [0JT8|j2|. In that 
problem D n had the form of an eddy-diffusivity operator: 


Dr sive (R) = ^4 eR 2 HR)4eSu(R) , 


R 2 dR 


dR 


(7.10) 


where h(R) is the eddy diffusivity which scales with R 
as a power law R^ h . Can we use this to guess a form for 
D n (R) in the present case? On the face of it the answer 
is no. Our evaluation ( |6.3| ) indicated that if we have a 
differential operator it operates on S n + 1 rather than on 
S n . On the other hand, once we assume that the leading 
order evaluation cancels in D n , the next order is again 
in terms of S n as in the case of passive scalar. In fact, it 
is not difficult to see (Appendix A) that the topology of 
the weakly linked diagrams for D n is identical (after the 
cancellation of the leading order) to the topology of the 
leading contributions for D n in the case of the passive 
scalar. We thus guess that for the aim of evaluation of 
the coefficient we can write D n (R) with an eddy-viscosity 
similar to (17.101) in which h(R) is found by comparing 
(pul) and MM: 


Cn = 


3-C2 


-1+1 I' 


4n(4 - C 2 ) 
3(3 - C 2) 2 


(7.16) 


Not having a theoretical value of (2 we can take for it any 
of the following values: (i) The K41 value ( 2 = 2/3. (ii) 
The experimentally “ accep ted” value <f 2 = 0.70. (iii) We 
can assume that Eq. (7.14) is correct also for n = 2. This 
gives ( 2 = 8/11 ~ 0.727. These different choices lead to 
the numerical values of shown in Table 1 under the 
rows labelled “model (i), (ii) and (iii)”. 

In the same table we also show the result of three avail¬ 
able experiments. It should be noted that all these ex¬ 
periments refer to Reynolds number Re a < 10 3 , and thus 
as discussed above, our asymptotic predictions may not 
apply. Corrections to our asymptotic behaviour are still 
too large to be ignored. This may be the reason for the 
deviations between the theoretical predictions and the 
experimental numbers. 


8. CONCLUDING REMARKS FOR PAPERS I-III 
AND THE ROAD AHEAD 

We believe that the theory that was presented in pa¬ 
pers I-III contains novel elements that are likely to re¬ 
main as cornerstones in the theory of the fine structure of 
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hydrodynamic turbulence. Since the approach is highly 
technical, we attempt in these concluding remarks to 
summarize first what are the essential elements of the 
analysis both from the point of view of technique and of 
physics. 

One fundamental issue that separates our approach 
from the traditional K41 approach is the explicit appear¬ 
ance of anomalous exponents with the ultraviolet cutoff 
length as the renormalization scale. It is intuitively clear 
that gradient fields cannot be insensitive to 77 since by def¬ 
inition a gradient measures differences on scales smaller 
than r], and this length may appear when one attempts to 
calculate the correlation functions of gadient fields. In¬ 
deed, an important step in our theory is the explicit iden¬ 
tification of the ultraviolet divergences which appear in 
the diagrammatic theory of correlations of gradient fields, 
and the computation of the largest anomalous exponent 
associated with such divergences. With this in hand we 
could go further in nonperturbative analysis and finally 
discover the mechanism for anomalous scaling. It is not 
completely evident that our route is either economical or 
unique. This route reflects the nature of our understand¬ 
ing and the path followed by us, and it is possible that 
there is a shorter path to the same results. Only future 
research will tell us about that. 

Let us summarize the technical aspects of our analy¬ 
sis. The foundation of the theory is the Navier Stokes 
equations. Starting from these equations we use a com¬ 
bination of renormalized (but order by order) analysis of 
diagrammatic perturbation series with some exact, non¬ 
perturbative considerations. We found that only the lat¬ 
ter allowed us to derive multiscaling. Order by order 
analysis could not take us out of the traditional K41 scal¬ 
ing laws. 

The scale invariant formulation of diagrammatic per¬ 
turbation theory for fluid mecahnics to all orders becomes 
technically tractable due to the Belinicher-L’vov trans¬ 
formation of the Eulerian velocity field. Without it all 
known approaches are either plagued by infrared diver¬ 
gences or are limited to low order perturbation terms. In 
paper I we developed a Wyld type diagrammatic tech¬ 
nique for BL velocity differences in r, f-representation. 
This presentation allows one to consider a p-tli order di¬ 
agram for the structure function S n (R ) as a set of “el¬ 
ementary” interactions involving p intermediate points 
with integration over their space-time coordinates. One 
important discovery of paper I is the property of locality 
which means that the major contribution to these inte¬ 
grals originates from a ball of locality which is a sphere 
of radius of the order of R. The second important fact 
which has implications on the structure of the perturba¬ 
tion theory is the property discovered in paper II which 
we called rigidity of diagrams. To understand this prop¬ 
erty one may use an analogy between diagrams and a 
mechanical system in which the propagators are substi¬ 
tuted by springs which are joined at the vertices. The 
point is that the dominant contribution to the diagrams 
for any object originates from the region of “mechanical 


equilibrium” of the complementary mechanical sytem. In 
the mechanical analog the external end of springs (prop¬ 
agators) are fixed at the external coordinates of the ob¬ 
ject discussed. The property of rigidity is very helpful 
in analysing the asymptotic behavior of many point ob¬ 
jects when some coordinates are much larger (or smaller) 
than others. It turned out that such asymptotic con¬ 
siderations for many point correlation functions play an 
important role in understandng the mechanism that leads 
to anomalous exponents. 

The main result of Paper II is the demonstration of the 
existence and origin of anomalous exponents in the an¬ 
alytical theory of turbulence.As in the theory of second 
order phase transitions there is a logarithmic divergence 
in the UV regime of ladder diagrams for correlation func¬ 
tions of velocity gradients. Anomalous behavior appears 
as a non-perturbative solution of some formally exact di¬ 
agrammatic equation which one obtains after resumming 
the ladder diagrams. Physically, anomalous behavior of 
correlation functions with two (or more) very different 
separation distances r«H is the result of a summation 
over a large number Jf = ( R/r) A of equally important 
channels of interaction of turbulent motions having scales 
between r and R. In paper II we succeeded to compute 
the anomalous scaling exponent A and found that A ex¬ 
actly equals its critical value A cr = 2 — £ 2 - This fact is of 
crucial importance for the nature of anomalous scaling 
in hydrodynamic turbulence. In the hypothetical case 
A < A cr one expects K41 scaling of the structure func¬ 
tions in the limit Re—> 00 . There can be only subcritical 
corrections to this, and such corrections have the viscous 
scale r] as the renormalization scale. In reality A = A cr 
and this opens up the possibility for the destruction of 
K41 for all values of Re. 

The analysis of the possible mechanisms for anomalous 
(non-K41) scaling of the structure function is the main 
topic of Paper III. This analysis is based on the exact and 
nonperturbative balance equations D n (R) = J„.(R) which 
follow from the equation of motion for the structure func¬ 
tions S n (R). These equations are a direct consequence 
of the Navier-Stokes equations in the statistical station¬ 
ary state in which dS n (R,t)/dt = 0. We succeeded to 
compute (in the limit Re—> 00 ) the viscous term J n (R) 
exactly: J n {R) = nJ2S n (R)/2S2(R). This result is the 
leading term among a few terms, and is the contribu¬ 
tion of strongly linked diagrams. It is proportional to 

RCn~(2 

as a direct consequence of the criticality of the 
theory in the sense that A = 2 — £ 2 - This contribu¬ 
tion is much larger than the contribution of non-linked 
diagrams, which are the result of Gaussian decomposi¬ 
tions. The largest of these next order terms is propor¬ 
tional to R^ n - 2 . In the hypothetical subcritical case when 
A < 2 — (,2 the latter term dominates the strongly linked 
terms and classical K41 scaling is then unavoidable. 

Because of the importance of this point we reiterate: 
because A attains its critical value 2 — £2 the dominant 
contribution to the viscous term J n (R) is proportional to 
S n (R)/S 2 (R) and not to S n - 2 (R)- This is one of three 
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findings which open the way to multiscaling in the the¬ 
ory of hydrodynamical turbulence. The second one is 
less solid: it is the assumption that the leading (stronly 
linked) contributions to the interaction term D n (R) can¬ 
cel exactly. We have only some arguments as to why 
this may happen. The detailed analysis of such a pos¬ 
sibility and the clarification of the possible relation of 
the cancellation with (hidden) symmetries of fluid me¬ 
chanics are important aspects of future research. What 
we propose now is that this cancellation is the only vis¬ 
ible mechanism to induce multiscaling in the theory of 
the structure functions. It is very possible that this is 
the only mechanism for multiscaling, and if we accept 
the opinion of the majority of the workers in the field 
of turbulence that experiments indicate the existence of 
multiscaling, we are forced to accept that there is a can¬ 
cellation of the strongly linked contribution to the inter¬ 
action term D n (R). Lastly, the third finding leading to 
multiscaling is the understanding of the important role 
of weakly linked contributions to S n (R) which scale like 
ftCn-i+Qj-Qi . j n a multiscaling situation this term dom¬ 
inates over the leading irreducible contribution which is 
proportional to R ^ n ~ 2+l ’ 2 . If one misses the weakly linked 
terms, K41 scaling reappears again. 

Besides issues of principle, there remains the actual 
calculation of the scaling exponents £ n . To accomplish 
such a calculation we need to know the coefficient in the 
evaluation of D n . From the “eddy viscosity” approxi¬ 
mation for D n we were led to the exponents shown in 
Table 1. We cannot say that the agreement with the ex¬ 
perimental data is amazing. The discrepancies may stem 
from various sources. Firstly, it is possible that the “eddy 
viscosity” approximation misses something important in 
the n dependence of the coefficient of D n . This is possi¬ 
ble, but may be not the main reason for the discrepancy. 
After all, in our theoretical developments we used very 
strongly the limit Re—> oo. This limit was assumed for 
example in neglecting the disconnected contributions to 
D n in favor of the weakly linked contributions. The dif¬ 
ference in the R dependence between these two types of 
terms (in the case of £>4 for example) is of the order of 
R 015 . We thus need at least 3-4 decades of inertial range 
to justify this procedure. It is very possible that the ex¬ 
perimental results are still suffering from the effects of 
subleading contributions, and only future analysis that 
takes such contributions into account may shed further 
light on the issue. 

One particularly pressing subject for future research 
is the asymptotic behaviour of n-order correlation func¬ 
tions. We found in the present paper that when two of the 
coordinates of F n { ri, r 2 ,... r n ) (say rq and r 2 ) coalesce, 
the correlation function is proportional to S 2 (rq). This 
is just one example of the asymptotic properties that are 
summarized under the term “fusion rules”. We need to 
understand the deep structure of the theory that is re¬ 
sponsible for this fusion rule, but that will also allow us 
to predict what happens, say, when three or more points 
coalesce. We guess that the n-point correlation function 


in that case will be proportional to F3, etc. Formally 
we need to develop the operator algebra that will auto¬ 
matically furnish all the needed fusion rules, and will be 
also compatible with multiscaling. We plan to propose 
elements of such a theory in a forthcoming publication 
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Finally, we comment on the relation and differences be¬ 
tween scaling in turbulence and scaling in better under¬ 
stood subjects like critical phenomena. Because of the su¬ 
perficial similarities (many body problems with strong in¬ 
teraction and scale invariance) there were many attepmts 
to apply formal schemes in the wake of critical phenom¬ 
ena to understand turbulence: renormalization groups, 
e-expansion, 1/d-expansion, 1/fV-expansions, and what 
not. If the approach taken in this series of papers turns 
out to be correct, this will mean that the theory of tur¬ 
bulence is significantly different from critical phenom¬ 
ena. Some elements reappear: sums of ladder diagrams 
contribute anomalous exponents, fusion rules are needed, 
and the smell of operator algebra is there. 

However that are at least two major differences: there 
exists flux equilibrium instead of thermodynamic equi¬ 
librium, and the interaction in the theory of turbulence 
is highly non-local because of the effects of pressure. In 
contrast, in critical phenomena it is sufficient to have lo¬ 
cal interactions that build up to global criticality because 
of the cancellation of energetic and entropic contribu¬ 
tions to the free energy. In turbulence, notwithstanding 
the nonlocality of the interaction it turns out that the 
diagrammatic theory in BL variables is finite order by 
order. In contrast, the perturbative analysis of critical 
phenomena leads to divergences that result (after renor¬ 
malization) in anomalous scaling. Thus the mechanism 
for anomalous scaling in turbulence must be different. 

Due to the flux equilibrium there is a global connection 
between the largest and smallest scales in the problem. 
A deep consequence of the flux equilibrium is the 2-point 
fusion rule that was discussed above as one of the corner¬ 
stones for multiscaling. In addition, flux equilibrium and 
the need to satisfy boundary conditions at the two ends 
of the energy cascade introduces the possibility of hav¬ 
ing the outer scale as the renormalization length without 
infrared divergences in order by order expansions for the 
structure functions. 

As a result of all these difference we do not have the 
fixed point structure with a small number of unstable 
directions that is so typical of critical phenomena. In 
some sense, the independence of the perturbative terms 
from a typical scale means that we have infinite number 
of marginal operators. The resummation of the pertur¬ 
bative theory results in a possibility of dressing these 
marginal operators, and there can be infinitely many in¬ 
dependent exponents. Whether such a concept can be 
turned into a computational scheme is a question for the 
future. 
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9. APPENDIX: WEAKLY LINKED 
CONTRIBUTIONS TO THE INTERACTION 
TERM 


In this appendix we discuss the weakly link ed co ntri¬ 
butio ns to the interaction term D n . In Eqs. (2.36) and 
(®> the integrals depend on the n + 1 -order correla¬ 
tion function T n +i (r,r,R) in which the two first coor¬ 
dinates r are special (they are either iq or R + iq and 
there is a derivative with respect to iq). In its turn every 
weakly linked contribution to F n+ i(ri, r 2 ,... r„+i) (or 
to T n+1 (r, r, R) = F n+1 (r, r, R, R, R) has two weakly 
linked legs (denoted in Fig.10a as x a ,Xb), connected to 
the body of the diagram via a one-propagator bridge. 
There are C% +1 = n(n + l)/2 weakly linked contributions 
to F n+ 1 in which the role of weakly linked legs is played 
by each pair taken from the n + 1 total number of legs. 
The C% +1 contributions to D n can be subdivided into 
three groups. The first group consists of just one term 
in which two special coordinates in T n+ \ are exactly the 
coordinates of the two weakly linked legs (r = r a = r^). 
The second group of terms in D n has 2 (n — 1) terms in 
which just one of the special coordinates in T n+1 (but 
only one of the two) is associated with a weakly linked 
leg. The second special coordinate is free to be asso¬ 
ciated with any of the (n — 1 ) legs of the body of the 
weakly linked diagram for F n + 1 . This body is an n-point 
object (see Fig. 10a) in which one leg is used to create 
the bridge. The last (third) group of terms in D n has 
= (n — l)(n — 2)/2 terms in which two special T- 
coordinates may be chosen from the coordinates of any 
(n— 1 ) free legs in the body of the weakly linked diagram 
for F n+1 . 

The first two groups of terms correspond exactly to 
the topology of the digrammatic representation of the 
interaction term D n in the problem of turbulent advec- 
tion of a passive scalar field T(r,f), and cf. section VB 
in p2| and section IIB2 in JTsj. Consider Fig.10 of p2[ . 
The dashed lines in this figure represent 2-point velocity 
correlators, and these are replaced in our case by wavy 
correlator lines. The wavy lines in the passive scalar case 
represent two point scalar correlators, and they remain 
as wavy lines in the present case. The fragment of the 
diagram in this figure which is placed to the right of the 
legs denoted by q 2 , 93 , <74 and to the right of the vertex 
between k and q\ must be interpreted now as a contribu¬ 
tion to the strongly-linked four-point velocity correlator. 
The last serves a a four- point body of a weakly linked 
contribution to a 5-point velocity correlator F$. We thus 
conclude that the topology of the diagrams for D n in 
the case of turbulent passive advection and the first two 
groups of weakly linked diagrams for D n in the case of 
Navier-Stokes turbulence is the same. 

It appears that the third group of C '^_ 1 terms which 
we described above forms a major difference between the 
problems of turbulent advection and Navier Stokes tur¬ 
bulence. In fact this group does not contribute. In the 


case of passive scalar this group is absent because of 
the zero value of the {Tv} correlator. It is remarkable 
that in the present case of Navier-Stokes turbulence this 
group cancels under the same condition of cancellation of 
the leading (strongly linked) contributions to D n . These 
terms may be considered (after severing the bridge to the 
weakly linked fragment) as strongly linked contributions 
to They must cancel if the scenario leading to 

multiscaling is assumed. 

The conclusion of this appendix is far from being triv¬ 
ial, and in some sense is very surprising. It says that al¬ 
though the passive advection problem is linear and local, 
whereas Navier Stokes in nonlinear and non-local (pres¬ 
sure!), it appears that if multiscaling is expected, the 
topology of the diagrams for D n is very similar in the 
two cases. If this is correct, it must be related to some 
deep symmetry. If so, the Eddy Viscosity Approximation 
used in section 7 may contain some essential aspects of 
the truth. 


* L’vov’s e-mail: fnlvov@wis.weizmann.ac.il. 
t Procaccia’s e-mail: cfprocac@weizmann.weizmann.ac.il. 
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FIG. 1. Diagrammatic representation of the 4-point 
Green’s function C/2,2 defined by Eq. (3.3). Panel (a): C/2,2 as 
a sum of weakly linked contributions G 2 ,2 (shown in Panels 
(b-d)) and strongly linked contribution C/|' 2 (shown on Fig. 
2). Panel (b):C/^2 is presented in terms of the three-point 
Green’s function C/2,1, defined by Eq. (3.1), and the dressed 
vertex A. Panels (c,d) give the diagrammatic representation of 
C/2,1 in terms of the propagators G and F 2 and the dressed ver¬ 
tices A and B. For the diagrammatic expansion of the dressed 
vertices A, B, and C - see Fig. 7 of paper II. Note the new 
notation of an empty little circle in the object V. This circle 
designates a vertex that can be either bare or dressed. 


FIG. 2. The diagrammatic presentation of the strongly 
linked contribution Qf 2 to the four-point Green’s function 
C/2,2, (see Fig.l). Diagrams 1 and 2 are the Gaussian decom¬ 
position of C/|* 2 , diagrams 3 and 5 are ladders with one rung 
(of two types), and diagrams 3 and 6 are the ladders with 
two and three rungs respectively. In contrast to the expan¬ 
sion in Fig.9 of paper II for G2, here one has rungs of three 
types. These are £3,1 like in diagram 5, E2,2 as in diagrams 
3 and 4, and rung Ei ,3 as the first rung in diagram 6. The 
diagrammatic expansion of £2,2 is shown in Fig.3, and the 
expansion of E3 ,i and Si,3 in Fig.4. The exact resummation 
of this series is shown in Fig. 5. 


FIG. 3. Diagrammatic expansion for the mass operator 
E2,2 which is found in the expansion of Q 2 ,2 shown in Fig. 
2. This mass operator is an important element in the system 
of equations for Q 2 2 and Qfi shown in Fig. 5. The principal 
cross sections are shown as dashed lines. 


FIG. 4. Diagrammatic expansions of the mass operator 
E3 ,i (panel a) and Si,3 (panel b) which appear in the ex¬ 
pansion of Q 2*2 shown in Fig. 2. These mass operators are 
the essential elements in the equations shown in Fig. 5. 


FIG. 5. The exact system of equations for the strongly 
linked contributions to G 2 ,2 (panel a) and C/3,1 (panel b) in 
terms of the propagators G and F2 and the mass operators 
E3 ,i, E 2,2 and Si,3. The first terms in the expansion of E m , n 
are shown in Figs.3, 4. This system results from the exact 
ressumation of ladder diagrams. 


FIG. 6. Diagrammatic representation of the weakly linked 
contributions to the four-point correlation function F4. Panel 
a: a weakly linked contribution to Fi(ro|;ri, * 2 , £ 3 , * 4 ) in 
which a one-propagator bridge is placed between (1,2) pair 
and the (3,4) pairs of legs. One has two more similar 
weakly-linked contributions with the bridge connecting the 
pairs (1,3) with (2,4) and (1,4) with (2,3). On the left of the 
diagrams 1 and 2 one find the objects G2 ,i and F3 resepec- 
tively. These objects are presented diagrammatically in Fig. 
lc and panel b of the present figure, respectively. On the right 
side of diagrams 1 and 2 one has the three-point objects that 
were designated as A and V. These objects can be found in 
Fig. Id and in panel b of the present figure. 


FIG. 7. Panel a: diagrammatic representation of the ir¬ 
reducible four-point correlation function as a sum of weakly 
and strongly linked parts. Diagram 1 is the weakly linked 
part FT 1 and diagram 2 is the strongly linked part FJ 1 . F" 1 
is presented in Fig. 6. Panel b: Exact presentation of FJ 1 . 
The elements appearing in this presentation are the four-point 
Green’s functions f/2,2 and f/3,1 on the left side and on the 
right side of the diagrams, and three different new four-point 
objects. The notation used to distinguish the three new ob¬ 
jects is with zero, one or two horizontal lines inside. The 
diagrammatic representation of the new objects is shown in 
Fig. 8. 


FIG. 8. Diagramamtic representation for the three different 
central blocks which appeared in Fig. 7. These are infinite 
order expansions in terms of contributions with increasing 
numbers of correlators in their principal cross section (i.e. 2, 
3, 4 and more. The four- point obejects in diagram 1 in each 
panel are the mass operators E2,2 and E3 ,i which are the 
dressed rungs of the ladders whose presentation is shown in 
Figs. 3 and 4. 


FIG. 9. An example of a fragment contribution to the right 
hand part of diagram 1 in Fig.7. 


FIG. 10. Diagrammatic representation of n-point correla¬ 
tion functions. Panel a: weakly linked contributions. Dia¬ 
grams 1 and 2 are the generalization of diagrams 1 and 2 for 
F4 shown in Fig.6 panel a. Diagrams 3 is an example of a 
weakly linked contribution with two one-propagator bridges. 
Such contributions do not appear in the case of F4. Panel 
b: Strongly linked contributions to F„, considering the legs 
designated as 1 and 2 special. In diagram 1 one has F4 on 
the left; in diagrams 2 and 3 one has G3 ,i and G2,2, the same 
objects that appeared in panel b of Fig. 7 for F4. 
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FIG. 11. Diagrammatic representation of the weakly linked 
contribution for F n shown as diagram 2 in Fig. 10, panel a. 
For the 3-point correlator F 3 (ro\x a ,Xb,x c ,Xd) we used the 
representation of Fig. 6, panel c and placed it here in the 
dashed ellipse which is denoted “left”. On the right we have 
an m-point object T> m which is a generalization for the case 
m > 3 of T >3 which appeared in Fig.l panel d. 
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